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Abstract 

In [TS], the author proved a number of multivariate elliptic hypergeometric integrals. The purpose of 
the present note is to explore more carefully the various limiting cases (hyperbolic, trigonometric, rational, 
and classical) that exist. In particular, we show (using some new estimates of generalized gamma functions) 
that the hyperbolic integrals (previously treated as purely formal limits) are indeed limiting cases. We also 
obtain a number of new trigonometric (g-hypergeometric) integral identities as limits from the elliptic level. 
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1 Introduction 

In [IB], the author proved and extended a pair of multivariate elliptic hypergeometric integrals conjectured by 
van Diejen and Spiridonov [4j, elliptic analogues of integrals due to Gustafson 9 . In a follow-up paper [5], 
van Diejen and Spiridonov proved "hyperbolic" degenerations of these integrals (after the hyperbolic gamma 
function of Ruijsensaars [18], and results of Stokman [22] at the univariate level). Unfortunately, the asymptotic 
estimates they had available were insufficient to derive these hyperbolic integrals as limits of the elliptic integrals; 
instead, they were forced to degenerate the elliptic proofs. The objective of the present paper is to put these 
(and other) degenerations on a sounder footing by showing that they are indeed limiting cases of the elliptic 
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integral. It is to be hoped that a better understanding of the relation between the elliptic and other integrals 
will lead to new results at various levels, both by clarifying how arguments at low levels extend to higher levels 
and by providing new results at low levels as limits of the elliptic identities. The latter hope has already been 
fulfilled to some extent; see, for instance, Corollary 15.51 below. 

A full study of degenerations of the elliptic hypergeometric integrals is beyond the scope of the present 
paper. Here, we focus only on the top level of each of the five types of limit (hyperbolic, trigonometric, elliptic, 
rational, and classical); degenerations within each type will be deferred to future work with F. van de Bult. 
We also do not consider in any detail the degenerations between types; rather, we obtain each case directly as 
a limit from the elliptic level. (Note, however, that in several cases our estimates are sufficiently uniform that 
one could view certain between-level limits as special cases of the limits from the elliptic level.) Similarly, we do 
not consider the discrete degenerations, a.k.a. hypergeometric sums; aside from finite sums (which have already 
been considered), the only interesting discrete degenerations appear to arise via lower-level (e.g., hyperbolic 
or trigonometric) integrals. Finally, we comment that at the univariate level, there are important relations 
between the trigonometric and hyperbolic integrals ( 22J, also 3 ); presumably, these have analogues at the 
multivariate level, but the existing techniques do not appear to apply. (However, see [21] for recent progress on 
a slightly degenerate, but multivariate case.) 

The plan of the paper is as follows. In Section 2, wc define the relevant gamma functions (elliptic, hyperbolic, 
trigonometric, and rational), and prove a number of asymptotic results relating them. In order to make the 
derivation of limit integrals as simple as possible, we have tried to make these estimates as uniform as possible, 
and the relevant error estimates similarly strong; since our techniques also apply equally well to higher-order 
multiple gamma functions, we prove our estimates in that level of generality. Even for the ordinary (i.e., r = 2) 
elliptic and hyperbolic gamma functions, our results are new, as the results in the literature were nonuniform 
(compare, for instance, Theorem 12.81 to its special cases Corollaries 12.51 and !2.9[ which were proved (for r = 2, 
and not uniform) in |18j). Corollary 12. 191 a uniform version of the q — > 1 asymptotics of g-symbols, may also 
be of independent interest. 

Section 3 considers the inequalities required to pinpoint where the integrands of interest are maximized. 
It turns out that these all follow from a single master inequality (Lemma I3.2[) , which in turn follows from 
an asymptotic analysis of an elliptic analogue of the Cauchy determinant. The result is an inequality stating 
that certain combinations of elliptic gamma functions are exponentially small unless their arguments alternate 
around the unit circle. It may be worth investigating other elliptic determinant and pfaffian identities to see 
whether they give rise to interesting inequalities. 

Section 4 begins the study of limit integrals with the hyperbolic case. This case turns out to be relatively 
straightforward, given the estimates and inequalities already established; in each case, a standard tail-exchange 
argument gives the desired limit. In fact, not only do we obtain the hyperbolic integrals as limits of the 
corresponding elliptic integrals, but moreover obtain exponentially small error estimates. We also discuss the 
extra arguments required to include the case in which the integrands are multiplied by appropriate abelian 
functions (the interpolation and biorthogonal functions of |16|). 

Section 5 considers the trigonometric (i.e., basic hypergeometric) case. It turns out that in addition to the 
integrals of Gustafson that originally motivated van Diejen and Spiridonov, there are additional limiting cases. 
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That these should exist is suggested by the fact that in the transformations of [16], only one side typically has 
a straightforward limit as p — > 0. It turns out, however, that if one breaks the symmetry of the integrand in 
a suitable way, one can arrange for both sides of the transformations to have reasonable limits. As a special 
case (violating our rule of considering only the top level of each case) , we find that not only are the ordinary 
Macdonald polynomials limits of the biorthogonal abelian functions, but in such a way that their orthogonality 
follows as well; in particular, the Macdonald "conjectures" (proved in [TT]) for ordinary Macdonald polynomials 
are limiting cases of the corresponding identities for biorthogonal abelian functions. 

The remaining two cases correspond to the case q — > 1 with p fixed. In the rational case, considered in 
Section 6, the parameters all tend to 1, and one obtains a hybrid of the hyperbolic and trigonometric cases. 
In particular, to obtain the full spectrum of possibilities, one must first break the symmetry of the integrand, 
and then do a tail-exchange asymptotic argument. The arguments are still fairly straightforward, although the 
resulting estimates are much weaker. 

Finally, in the classical case, considered in Section 7, the parameters behave in such a way that the integrand 
is exponentially small unless certain inequalities are satisfied. This gives rise to multivariate analogues of the 
ordinary beta integral; notably those of Dixon [5] and Selberg 19J. Curiously, the most natural forms of these 
limits remain elliptic in nature, with the integrand involving powers of theta functions, although these can be 
removed with a suitable change of variables. 

The author would like to thank P. Forrester, J. Stokman and F. van de Bult for motivating conversations 
regarding the trigonometric and hyperbolic cases, and R. Askey for suggesting the use of the modular transfor- 
mation to derive classical limits (as in [13]), which led the author to consider the paper [12]; the author would 
also like to thank an anonymous referee for pointing out that the original version of Theorem 14.71 was badly 
stated. The author was supported in part by NSF Grant No. DMS-0401387. 

Conventions 

We use standard conventions for big O notation in uniform estimates; that is, we say that f(v, z) — 0(g(v, z)) 
as v — > + uniformly over the region z S D(v) (where D(v) is a set depending on v) if there exist constants 
S > and C > such that 

\f(v,z)\<C\g(v,z)\ 

whenever < v < 8, z 6 D(v); and similarly for limits v — > ±oo (which can be viewed as limits ±l/v — > + ). 

All logarithms are taken on the principal branch, with branch cut along the negative real axis. All powers 
are determined correspondingly; in particular, square roots are chosen with positive real part. (The obvious 
exception is which is always taken to have positive imaginary part; we refrain from denoting it by i so as 

to retain i as an indexing variable.) 

Finally, for a real number x, {x} denotes the fractional part of x; i.e., the unique representative in [0, 1) of 
x + Z. 
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2 Asymptotics of multiple gamma functions 

For integers r > 0, let (w; x; u>% , . . . , u> r ) denote the function 

} (w;x;gji, ...,u) T ) 



(r) x . e ( xw ) 



rii< 4 <r( e (^ w ) - !)' 

where e(x) := exp(27r^/— Tx), and define a family of polynomials Pn (x;u>i, . . . ,co r ) in x by the Laurent series 
expansion 

(r) (w;a;;wi, . . . ,w r ) = w n P^\x] uo ll . . . , w r ) 

n> — r 

valid in a punctured neighborhood of w — 0; we also set 

Q« (*;«!,... ,w r ) :=P^+ 1 )( ; r;a; 1 ,... ) c i ; r) -l). 

By convention, if n is omitted, then n = 0. Note that these are related to the multiple Bernoulli polynomials 

of [T5] by 

(2tt\/— T) n 

■fn (^J^l) •••,W r ) = 7 j ri-B r , n +r(a;|Wl, . . . ,W r ). 

(n + r)\ 



Note the special cases 



P(°)(i;) = l 



2 , ,2 , ,2 
2 



p(a)( j = 12(,-^/2-c 2 /2)^-^-c 



and 



)(°>(x; ) = -(* + 1/2) 

-12(x - oj/2 + 1/2) 2 + w 2 + 1 



24w 

\2 , ,.,2 , ,.,2 



„ (2) , . (x-wi/2-w 2 /2 + l/2)(-4(x-w 1 /2-w 2 /2 + l/2) 2 +w 2 +w| + l) 
Q ^^) = ■ 

We will also need a third polynomial 

R (r \x;uj 1: ...,w r ) = Q (r) (x;a;i, . . . , w r ) + ip (r) (x; wi, . . . ,u r ) 

Q { - r \x; wi, . . . , w r ) + (3 (r) (x - 1; wi, . . . ,w r ) 



and note 



i? (0) (x;) =-x 
^)(,;.)=- 12( — /2)2+ " 2 " 2 



(2) (x-w 1 /2-c; 2 /2)(-4(a;-c l ; 1 /2-w 2 /2) 2 + u; 2 +w 2 -2) 
(^^i'^ 2 ) = 24^ ' 



We can then define the hyperbolic gamma function vj^ as follows. For all integers r > 1, and < < 

El<,<r ^(^i)' 

rl r \x;u>i,u 2 , ■ ■ ■ ,u r ) := exp(PvJ </> (r) (w; x; u%, . . . ,w r )-^), 

where by the principal value integral PV/ R we mean the average of the integral over two contours agreeing with 
R away from 0; one contour passes to the left of 0, and the other passes to the right of 0. This differs slightly 
from the corresponding definitions of hyperbolic gamma functions in the literature; in particular, we have 

Fj l r \x;uj 1 , ...,u r ) = S r (x\ui, . . .,w r ) ( ~ 1) 

where S r is the multiple sine function (see, for instance, [12j). and 

/ 2\ 

l£ '(x;ui,u} 2 ) =G(-cjiV-1,-^2V-1;x-cji/2-cj 2 /2) 
in terms of Ruijsenaars' hyperbolic gamma function [18] . Since 



PV 

and 



( r \w,x]L>i, . . . ,u r )— = ^Pv/ (^ r \w;x;u>i, . . . ,w r ) - <ft r \-w;x;uji, . . . ,u r j) — 



PVy w k dw = 

for integers k < — 1, we can express this as the ordinary integral over R of 
e(xw) + ("l) r_1 e((Ei<i<r^ ~ x ) w ) 



2^ TTi — 1) 

lll<i<rV V J J l<fc<(r+l)/2 



-2k 



or by symmetry as twice the integral over [0, oo), thus recovering the definitions of [18] and |12j . 
When r = 1, we have 

e(wx) dw ,„ . , ... 

Jr (e(iuw) - 1) w 

for < 3(x) < 3(o;) (the branch with value log(2) at x = w/2), which then gives an analytic continuation of 



1^ to all x, uj such that w^fl: 



r^ 1 \x;u>) — 2 sin(7rx/cj). 



(This is a fairly standard definite integral; it can be shown, for instance , by moving the contour infinitely far 
to the left (for x/oj) > 0; if x/uS) < 0, the contour should be moved to the right, and the case —x/w 
real follows by analytic continuation), and observing that the resulting sum of residues (compare Theorem [ 
below) is 

ny/—l(x/Lj) — 7rV— 1/2 — e(—kx/u))/k = log(2 sin(7ra/u;)), 

fc>i 

as required.) In general, we have 

T^ r \x + w r ] Ui, ...,w r ) = tI[\x]uj u . . . ,cj r )r^ r_1) (x; uj 1i . . . ,uv_i) 
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(r) 

in the domain of definition, and thus by induction have a meromorphic continuation to all x. Since was 

(r) 

defined via its logarithm, we in fact have a nearly canonical choice of branch for logl^ , more precisely, for 
each cuq with ^v(^o) > 0? there is a unique analytic continuation of 

f duu 
logr^ r) (a;; wi, ...,w r ) := PV7 (w; x; wi, . . . ,u r ) — 

Jr w 

to the set (uq; u)\, . . . , uj r ) obtained from C by removing the countable unions of rays 

w R>o + ^ ^iZ>i and w K<o + ^ WjZ< . 

l<i<r l<2<r 
(r) 

That is, for each zero and pole of 1^ , we cut along a ray in the direction ±wo, as appropriate. Also of 
importance is the analogous analytic continuation of 

log((a:/wb) ( - 1)r rM(a!;w 1 ,... ) w r )) 

to the domain C'^ r \u!o;u>i, ■ ■ ■ ,uj r ) which differs from C^^ujq; oj\, . . . ,u> r ) only in that the ray woR<o has not 
been cut; this continuation exists since the zero/pole at has been cancelled. We also by convention take 
rjf\x; ) = — 1, to make the functional equation valid for r = 1 as well. For the functional equation to hold for 
the logarithm, we must take 

logr^ 0) (x;) := -7rV^Tsgn(9(a;/cjo)), 



defined on the set C (0) (w ; ) = C \ lo q R. 
Remark. Using the fact that 

(cx; clji,..., cuj r ) = (x;u>!,...,cj r ) 

we can further extend rj^ to arbitrary u>i, so long as there exists a constant c with 3(co;i) > for all i. We 
will not be using this extension in the sequel, although the ability to rescale within the upper half-plane will be 
quite useful in the proofs. 

In [HI Prop. 5], Narukawa derived a product expansion for based on the observation that the integral 
over a suitably chosen sequence of contours = a can be made to tend to as a — > oo, and thus the integral 
expands as a sum of residues. We will need a more precise form of the bound on the integral. 

Lemma 2.1. Fix e > 0, and let a, Q(x), LOi range over the domain < < $>(u)i), and for 1 < i < r, 

^s(u>i) > and \a — n/u>i)\ > e, all integers n > 0. Then 



{r) dw 



<a~ l 



cxp(— 2nadt(x)) exp(— 2ira'tR.(x — J2i w *)) 



where C(x) = min t 2(a,z)>oi | e ( j/) — 1 1 . In particular, the integral is uniformly 0(exp(—2na'Si(x))) over any compact 
subset of the domain. 
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Proof. We have the estimates 

|e(ioa;)| = exp(-2vr3fj(w)3(x)) cxp(-27ra3f?(x)) 
\iu\- 1 = l/V^M 2 +a 2 < l/o 
|e(«xji) - I!" 1 < C(elwil)- 1 

\e(wui)- < C(e|w i |)- 1 exp(27r3fJ(w)3(w i ))cxp(27ra3?(w i )) 

where we note that |wcjj — n| > e|wj|. If we use the third estimate for 3?(w) > and the fourth estimate for 
Sft(io) < 0, we obtain the stated bound. □ 

Narukawa then uses the fact that for $l(x) > max(0, 3?(X)i w »))> tne above bound tends to as a — > oo; for 
our purposes, it is more convenient to fix a and obtain an asymptotic series. 

Theorem 2.2. Let a, arg(x) 7 Wj range over fAe domain a > 0, 9f(e(— arg(x)/27r)a>i) > /or 1 < i < r. Then 

as x — > oo, 

-7iV^IP (r) (a;;wi, . . . ,uy) - Res t0=y ((/> (r) (w;x; wi, . . . ,w r ) — ) + logi^ r) (x; u>i, ...,u> r ) 

0<9(e(arg(a;)/27r)j/)<a 

= 0(cxp(— 27ra|x|)), 

uniformly over any compact subset of the domain. 

Proof. Clearly replacing arg(:r) by arg(x — 6 ^ oji) for < b < 1 will have no effect on the validity of the bound. 
We may thus take 

x = e{r)\y\ +6^^. 

i 

Since 

logT^ (e(4>)\y\ + fc^Wijwi, . . . ,w r ) = \ogT%\\y\ + b e(-<£)wi; e(-<l>)<Ji, . . . ,e(-<t>)u r ), 

i i 

we find that uniformity in a, r, u>i will follow from uniformity in a, Ui with r = 0. 

Thus assume = k^i^K^) f° r < 6 < 1. For every point in the domain, there exists a' > a and 

e > such that the previous lemma applies, and such that there are no poles with imaginary part in (a, a']; 
by compactness, we can cover the domain by a finite number of such choices. Since (9(exp(— 2na'$t(x))) = 
(9(exp(— 2wa$l(x))), it suffices to consider the case a' = a. The result follows by residue calculus. □ 

Remark. Note that the cut lines for log vj^ can be taken along any direction lu in the convex cone generated 
by Ui,...,u n and the argument is valid for x in the complement of the cones 

u>i + ■ ■ ■ + ay + M<o (u>i , u) r ) an -d — K<o (o»i , . . . , u> r ) , 

and not approaching oo parallel to the boundary of the cone. 
In particular, we obtain the following estimate. 
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Corollary 2.3. Let a, a,rg(x), u>i range over the domain < a < mini 3(— e(arg(x)/27r)/cji). Then as x — > oo, 
we have the estimates 

-7r\/^TP (r) (x; wi, . . . , uj r ) + logl^O; ux, ...,co r ) = 0(exp(-27ra|a;|)), 
ny/^lP^i-x; u) X , ...,u r ) + ]ogT^\-x; ui, ...,w r ) = 0(exp(-27ra|a:|)), 

uniformly on compact subsets. 

Proof. The first estimate follows immediately from the theorem. The second estimate then follows using the 
facts 

pM (x;wi,...,w r ) = (-iypM £wj-a;wi w r ) 

i 

and 

logr^ r) (a;;cji, . . .,u) r ) = (-l) r_1 logT^ r) (^ uj, t - x; Ui, . . . , u r ) 

i 

□ 

Remark. For r = 2, this result is essentially due to Ruijsenaars |18[ App. 2]. In general, on the domain 
9f(— e(arg(a;)/27r)/a;i) > 0, it should be possible to improve the error term to 

0{\x\ r ~ x exp(-27r|a;| min3(-e(arg(x)/27r)/cj. i )), 

i 

by taking a = mim 3(— e(arg(x)/27r)/wi) and bounding the leading residues. In particular, if the poles with 
minimum imaginary part are all simple in the given compact subset of parameter space, then one easily has the 
bound 

0(exp(— 2n\x\ min3(— e(arg(x)/27r)/c<;j))) 

i 

on their residues, and thus on the error term. 

We will also need an extension of this to the case that one of the moduli tends to 0. For convenience, write 

lh\ x 'i w i> • • - ><*V) = -7rV^TpW(x;a;i, ...,w r ) + logr^ r) (x; oji, . . .,co r ). 

Theorem 2.4. Let a, arg(z), uji,a,(3 range over the domain 

0<a< min 3(— e(arg(x)/27r)/wj), Q ! (e(— arg(x)/27r)/w r ) > 0. 

l<?'<r— 1 

Then as x — > oo, v — ► + , 

7h\x + vaw r ',wi, ■ ■ .,w r -i,vu r ) - ^\x + v/3u! r ;oJi, . . .,w r -i,vu r ) - (a- f3)j^~ 1 \x;uj 1 , . . .,w r -i) 

= 0(v exp(— 27ra|a;|)), 

uniformly over compact subsets of the domain. 
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Proof. For r = 1, this is immediate, so take r > 2. As before, we may restrict our attention to the case 
i<«r-i Then the left-hand side can be expressed as 

(e{vau> r w) - e(vf3uj r w) \ (r _ iv .diu 

I3t(w)= a \ e(vuj r w) - 1 J w 

The desired estimate then follows as in Lemma I2TT1 using the fact that as v — * + , 

1 / e(vau) r w) — e(v(3uj r w) 

- 7 \ ; (a- P) 

v \ e(vu} r w) — 1 

is 0(exp(e|3?(ii;)|)) for any e > 0, uniformly in w and on compact subsets of parameter space. □ 

Corollary 2.5. Let x, . . ,u r) a, (3 range over the domain > 0, x G C^' r_1 ' ) (a; r ;cj 1 , . . . ,cj r _ 1 ). Then 

as v — > + , we have the estimate 

log -- (a-/3)logI^ ^x;^,...,^.!) =0(u), 

uniformly over compact subsets of the domain. 
Proof. Since 

C( r ~ i:) (w r ;wi, . . . ,uv-i) = C (r) (w,,; Wi, . . .,u r ), 

there is a canonical choice of branch for the left-hand side, and the proof of Theorem 12.41 gives the desired 
estimate on compact subsets of the complement of the cones 

uji ^ + w r + R<o(u;i, . . . , uv) and — K<o(o;i, . . . , u r ). 

(r) 

Since the estimate is consistent with the functional equation for 1^ , the result follows. □ 
Remark. In both cases, this can be extended to a uniform asymptotic series 

T ( £ ) {x + vau r \ui,...,uj r -.-i,vu r ) _ ^ B k (a) - B k ((3) ( d \ k ~ 1 r (r-i). 



log-^y— — -= ^ ^ logr ^ '(x;ui,---,^ 1 )+0(v ), 



with error 0(w™ exp(— 27ra|x|)) if also x — > oo; here B k {x) is the fcth (ordinary) Bernoulli polynomial. 

In fact, with care, we can give an estimate valid on the larger domain C'^ r ~^' {u r \u)x, . . . , w r _i), and thus 
in particular in a neighborhood of 0. The point is that we can identify the zeros and poles of T fl and 
that give rise to the cut line <x> r lR<o above, and using the ordinary gamma function, cancel them out. The 
resulting asymptotics can then still be computed via Stirling's formula. It will be convenient to use a slightly 
renormalized form of the ordinary gamma function; we define 

with the convention that u = 1 if omitted, and the usual convention on multiple arguments. As a justification 
for this convention, note that the reflection identity for the ordinary gamma function becomes 



T r (x, w — x;u) = Y^\x;uS) 1 . 



9 



We also let logr r (a:;a;) denote the standard branch on C \ wK<o, and note the following version of Stirling's 
formula. 

Theorem 2.6. Let x ^ M< . Then for all m > 1, 

lo g r r ( a O = (zlog(z)-z)-~log( a; )+ £ ^^- +O(d(x,R< ) 1 - 2m ), 

l<i<m ^ ' 

uniformly in x as d(x,M<o) — * oo. More generally, for a ranging over any compact subset o/C, 

logT r (x + a) = (xlog(x)-x)+B 1 (a)log(x) + V t^l^M + (d(x, R< )' m ), 

'-^ tit + l)x l 

l<i<m 



uniformly as d(x, M<o) — > oo. 



Proof. The claim for general a follows from the claim for a = by straightforward algebraic manipulation. For 
a = 0, we observe that by [121 §8.4] or J5J Thm. 1.4.2], the error term is 

B 2 m ({t})_ dt 



2m(2m - 1) 2m J (x + t) 2m 
The first term certainly has the correct asymptotics; for the second term, we have 

i r B ^({t}\ dt=0{ r lx+tl -2m dt y 



2m J (x + ty 

For 9ft(x) > 0, \x + t\-' 2m < (\x\ 2 + t 2 )-' m , and thus the integral has order \x\ 1 ~ 2m . For R(at) < 0, the integral 
is still bounded above by 

/oo 
(3(x) 2 +t 2 y m dx = 0(\^s(x)\ 1 - 2m ). 
-OO 

□ 

Remark. More generally, [TU Ex. 8.4.4] gives the error term 

(-l) m+1 B m+1 (a)x- m 1_ Z" 00 B m+1 ({t-a}) 

2m(2m+l) m+lj (x + t) m+1 

for < a < 1. 

Corollary 2.7. Let x,uj,a,f3 range over the domain x£C \ cjIR<o. Then as v — > + , 

l0 /r(* + va^) E (-l)'(^(a) - ^(^(.T/t^) 1 -' 

+ O(d0r,^M< )-"V"), 
uniformly in x and over compacta in w, a, (3. 
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Theorem 2.8. Let x, uj\,. . . ,uj r , a,(3 range over the domain $s(u)i) > 0, x € C'^ r 1 ^(w r ;wi, . . . , ui r -x). Then as 
v — > + , we have the estimate 

1^ (a; + vauj r ; uj\, . . . , u) r -i, vuj r ) T r (x + vau r ;vui r ) 

lOg — (— 1) lOg ; 

^'{x + v(3uj r - u u . . . , u r -i, vuj r ) V r {x + vj3u r ;vuj r ) 

-(a-p) {\agT£- x \x\wi, . . .,w r _i) - (-lflogfc/tK*.)) = 0(«), 

uniformly over compact subsets of the domain. 

Proof. We first observe that 

logl^ (a; + vau r ; u>x, . . . , w r -i) — (— l) r logr r (x + waw r ; !)w r ) 

and 

logT^^ix; u>i, . . .,u> r -i) - (-l) r log(a;/i;w r ) 

are analytic on the given domain, and thus the overall left-hand side is analytic. Moreover, the stated estimate 
holds on the smaller domain C^ r ~^\u) r ; lo\, . . . , ui r —x). Using the functional equation, we may immediately 
extend this to the full domain C'( r_1 )(u; r ; w\, . . . ,u; r _i), except when r — 2, where the point x = must still 
be excluded (since the only points related to x — via the functional equation are also in the cut set). But in 
that case, we may simply use Cauchy's theorem to deduce a uniform estimate on a neighborhood of from a 
uniform estimate on the boundary of the neighborhood. □ 

Remark. More generally one has an asymptotic series in which the coefficient of the kth term depends on the 
(k — l)-st derivative of 



If r = 2, we have 



logI^ r x \x;ui, . . .,u>r-i) - (-l) r \og(x/vu) r ). 



. \ t„... \ _ l0g(wi) - l0g(w 2 ) - l0g(27T) 



lim logT^ J (x; wi, w 2 ) - logr r (x; w 2 ) . 

x^O Z 

which gives rise to a nice corollary by taking x = (3 = above. 
Corollary 2.9. UBS As v -> 0+. 

(2) 

log]?^ ' (t)aw 2 ;wi,w 2 ) — logr r (a) — (a — 1/2) log^Triiu^/^i) = 0(v) 
uniformly over compact subsets of the region Si(a>2) > 0, a S C. 

Now, consider the elliptic gamma function, defined as 

r^(z; Pl ,p 2 , ... )Pr )= J] (1 - Pl 1 + 1 p k 2 2 + 1 ■ ■■P h r r + 1 /z)(l - P k M 2 ■ ■■P k r r 4- 1Y 

0<fcl,fc 2 ,...,fcr 

For \p\P2 ■• -Pr\ < M < 1, we have 



logri >(z;pi,p2, ..., Pr ) = ^ -= — — : 

i<- t K ILKikA 1 - ~ Pi ) 
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this, then, for any ui with 3(wo) > 0, defines a branch of logr e on the region Ce (wolPij • ■ ■ ,Pr) obtained 
from C* = C \ {0} by removing the countable union of logarithmic spirals 

Z>1 Z>1 Z>1 / ^ N j Z<Q Z<Q Z<Q / TTT) \ 

Pi" P2" ' ' 'Pr~ e(w R>o) and p x - p 2 - • • <p r - e(cj R<o)- 
(We define a region Ce analogously.) And, of course, we have the functional equation 

logT e (r) (p r z;pi,p 2 , . . . ,p r ) - logT e (r) (z;pi,p 2 , ...,p r ) = logr e (r_1) (z;p 1 ,p 2 , . . . ,p r -i), 

with 

lo g ri°)(z;) = log(-l/z), lo g ri°)(-l;)=0. 

Narukawa [T2l Theorem 14] gives the following "product" expansion of the elliptic gamma function in terms 
of the hyperbolic gamma function: 

-27r>/=TQ (r) (a:;wi ) . . . , w P ) + logrJ p >(e(x);e(wi), . . . , e(w r )) 

= -TiV^IpM (a: + fc + 1; wi, . . . , bj r ) + logics + k + 1; u u . . . , u r ) 

0<k 

+ TrV^lP {r} {x - k; uji , . . . , w r ) + log (x - k; oji , . . . , w r ) 

0<k 

Note that each term in the infinite sums converges uniformly exponentially to as k — > 00, so the sums converge 
uniformly and absolutely. 

Using this expansion and the asymptotics of the hyperbolic gamma function, we obtain the following esti- 
mates for the elliptic gamma function. First, the hyperbolic limit pi, . . . ,p r — ► 1. 

Proposition 2.10. Let A, . . . ,u) r , x range over the domain 

< A < min 9?(-l/wi) - \$l(x/u)i)\. 

l<i<r 

Then as v — > + , we have the estimate 

-2-K\/~^lR^ r \x; vuji, . . . , VLu r )+logT^ r \e(x);e(vu!i), . . . , e(vuj r ))-logT^ (x; vui, . . . ,vuj r ) = 0(exp(-27rA/v)), 
uniform over compact subsets of the domain. 

Remark. If we do not care about the choice of A, the constraint on the domain is simply 

|S>(x/«i)| < 3(-l/wi),l< i < r; 
this is a parallelogram, two of the vertices of which are ±1. 
Proposition 2.11. Let A, uj\, . . . ,uj r , x range over the domain 

0<A< min min(9f(— x/u>i), $s((x — l)/u>i)). 

l<i<r 

Then as v — > + , we have the estimate 

-27r\/^TQ ( '' ) (a;;uw 1 , . . . , vuo r ) + logr e (r) (e(a;); e(vw{), . . .,e(vw r )) = 0(exp(-27rA/u)), 
uniform over compact subsets of the domain. 
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Proof. Using the "product" expansion and the identity 

we can express each left-hand side as a sum over functions to which Corollary 12.31 applies, giving the desired 
uniform asymptotics. □ 

Using this limit, we can obtain the following bound. 

Proposition 2.12. Fix a compact subset S of the set of u>i such that < and constants < e, C\, C%. 

Then as v — > + , we have the following estimate, uniform over the product of S with the region —1 — vC\ < 
9?(x) < vC\, \^s(x)\ < vC2, apart from a hole of radius ev around every pole of the left-hand side: 

rP(e(x); e(wwi), . . . , e( Wr )) ±1 = 0(e(±Q< r >(x; vu u ..., «w r ))). 

Proof. We consider the + case; the — case is completely analogous. Choose 1/2 < D < 1 and < a < 
minj 9f(— l/u>i). Proposition ^. 101 gives 

e(-Q (r) (x; vwi, . . . , vu) r ))T^ (e(x); e(wi), . . .,e(w r )) = 0(e(ip (r) (a;; wi, . . . , w I .))r^ r) (a;; vwi, . . .,voj r )) 
away from the poles, in the subregion — D < §t(%) < vC\. Since 

e(^-P (r) (£; v^i, . . . , vu) r ))T^ (xyvwx, ■ ■ - ,vu r ) 

is uniformly bounded in that region (we have excluded neighborhoods of the poles, and it converges uniformly 
to 1 for — 5ft(x)/w large), we have the uniform estimate 

e(-Q (r \x; vuix,..., voj r ))T^ (e(x); e{vw x ), . . . , e(vcu r )) = 0(1) 

in this region. A similar application of Proposition ^. 101 gives 

e(-Q (r) (x + l;vu>i, . . . , voj r ))T^ (e(x); e(voJi), . . . ,e(vu r )) 

= 0{e{^P {r) (x + l;vwi, . . . , vuj r ))T ( h r} (x + l;vu x , . . ., vw r )) 

uniformly over the subset — vC% < $t(x) + 1 < D. Since 

Q^(x + 1; vuji, . . . , vu> r ) — Q( r '(x; vu>i, . . . , vuo r ) — P^ r \x + 1; vuj\, . . . , voj r ), 

we find 

e(-Q {r) {x]vuji, . . .,vu> r ))T^(e(x);e(vuJi), . . .,e(vu) r )) 

= 0(e(~P^{x + 1;vlj 1; VLU r ))T { h r) (x + 1; vljx, vu r )) = 0(1) 

on this region as well. □ 

(r) 

Similarly, we can obtain asymptotics of r e in the "rational" limit p r —> 1, p lv . .p r —\ fixed. 
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Theorem 2.13. Let x, p\,. . . ,p r -i, w r , ct, (3 range over the domain < \p\\, ■ ■ ■ , |pr-i| < 1, 3(w r ) > 0, and 

x e e- 1 (Ci r - v >( ! Jr',Pi, ■ ■ ■ ,Pr-i)) U KR> n e-^C'^-V^px, . . . ,Pr-l)). 

Then as v — > + , we have the estimate 

log ^ r \z{x + wa^r);Pi, ■ ■ ■ ,p r -i,e{vuj r )) _ , ^y^ r r (x + vau) r ;vu} r ) 
vt\e(x + v0LJ r y,p!, . . . ,p r -i,e(vLJ r )) r A x + vPu) r ;voj r ) 

-(a -(3) (loglf- 1 )^);?!, . . . , Pr ^) - (-l) r \og(x/vu; r )) = 0(v), 
uniformly over compact subsets of the domain. 

Proof. If we expand the elliptic gamma functions via the product representation and group corresponding terms, 
we find by Corollary |2 . 51 that all but one term of the result is uniformly 0(v). Moreover, it follows from Theorem 
12.41 that the coefficient of v in the estimates is exponentially small as k — > oo, and thus the error terms are 
summable. The only surviving term can be estimated using Theorem 12. 8[ giving the desired result. □ 

Corollary 2.14. As v ^ 0+. 

log ri 2) (e(at r );p,eM) = ^ + _ ^^.^/^ + Q{v 

uniformly over compact subsets of the region $s(ui) > 0, < \p\ < 1, a G C. 

Similarly, one has the following. 

Corollary 2.15. As v ^ 0+. 

Tt\e(avu J ) Z ;p 1 ,...,p r . 1 Mvu J )) = + 0(t>)) 

Te ' (e(f3vu:)z]pi 1 . . . ,p r -i, e{vuj)) 

uniformly over compact subsets of the domain < \p\\, ■ ■ ■ , 1| < 1; > 0, \pi ■ ■ ■ p r -i\ < Z < 

mini<i< r _i \pi\~ 1 , z ^ 1, a, (3 € C. 

Corollary 2.16. As v — ► + , f/ie function 

Te V \e{vauj)z]pi, ...,p r , e{vuj)) 
rt\e(v(3uj)z;pi, . . . ,p r ,e(vu>)) 

is uniformly bounded over compact subsets of the domain < |pi|, . . . , |,Pr-i| < lj > 0, |pi • • ■ j?r— 1| < z < 

1, ex, [3 £ C. On compact subsets of the domain < |pi|, . . . , |p r -i| < 1> > 0, \p\ < z < 1, a, (3 £ C, i< is 

uniformly 

O / w min(SR(a-^),0)-j_ 

There are corresponding estimates for the trigonometric gamma function 

0<fc v lHJ 

with the corresponding analytic continuation of its logarithm (the branch with logT t (0; q) = 0). 
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Lemma 2.17. Let z, oj, a, (3 range over the domain > 0, < \z\ < 1. Then as v —> 0, we have the 

estimate 

T t (e(vauj)z:e(vu!)) , „. , . , . . 

lo § f a ( - « - log 1 - z = OH, 

r t (e(^w)z; e(TO)J 

uniformly over compact subsets of the domain. 

Proof. Indeed, for \z\ < 1, the left-hand side is given by the sum 

Ez k ( e(kvaui) — e(kvf3ui) \ 
T { 1 - e(kvuj) + 01 ~ ' ) 

k>l v v ; 7 

and the quantity in parentheses is uniformly 0(v exp(efc)) for all e > 0. □ 

This gives rise to a trigonometric analogue of Theorem 12.131 

Theorem 2.18. Forr > 1, let z, p\,. . . ,p r -\, u>, a, (3 range over the domain < \pi\, . . . , \p r — 1| < 1, > 0, 

and z £ Ce ■ ■ ■ iPr-l)- Then as v — > + , we have the estimate 

log r e (r) (e(t;a^)2:;pi, . . . ,p r ._i,e(>^)) _ ^ r^(e(^a^)z; e(ro)) 
r e (r) (e(w/3o;)2;pi, . . . ,p r _i, e(w)) T t (e(v(3uj)z; e(w)) 

- (a - /3) (log it" 1 ) (2; p 1; . . . ,p r _ x ) - (-If log(l - *)) = O(w), 
uniformly over compact subsets of the domain. 

Proof. For \p1p2 ■ ■ 'Pr-il < z < rninj (pil^ 1 , this follows immediately from the expansion 

logr«(z;pi, . . . ,p r ^,q) = - J2 r *^ 1+1 ' ' -P k r-i +1 h^) + Mr 1 ^ 1 • ■■P h r - 1 1 z;q) 

0<fci ,&2 ,. . . ,/c r - 1 

and the asymptotics of T t . 

The general case follows from the functional equation. □ 



Comparing this to Theorem 12.131 gives the following result, a uniform version of the results of the appendices 

of pig. 

Corollary 2.19. Let x, uo, a, [3 range over the domain > 0, 

x € (C \ (Z + cjR>o)) U wM> . 

Then as v — > + . we /iawe i/ie estimate 

r t (e(x + vauj); e(vuj)) T r (x + vauj; to) , , 

lo § rTTTTl 1 — \\ ~ lo § r ? , 5 V ~ a ~ ^ ( lo S 1 ~ e W - logO/TO ) = 0(v), 

i t {e(x + vpuj);e{vuj)) L r (x + vpui; to) 

uniformly over compact subsets of the domain. 

Remark. Note that the validity of the theorem for r = 2, |p| < z < is enough to give the corollary in 

general, and in turn give the theorem in general, without using the functional equation. This also implies that 
Theorem 12.131 and its corollaries continue to hold even without the constraint < \pi\ on the domain, and 
further implies that Lemma |2 . 1 71 holds on the domain z ^ e(ajR>o). 
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Corollary 2.20. \10lj Let u, a range over the domain > 0, a G C. Then asi>~> + , we have the estimate 

log r( \ = T^— + (o - 1/2 log(27ro«/V=i) + O « , 

l r (a) 12w 

uniformly over compact subsets of the domain. 

3 Generalized triangle inequalities 

When taking limits of elliptic hypergeometric integrals, the first step is naturally to determine which part of 
the contour makes the most significant contribution to the integral. We first note the following consequence of 
Proposition 12.12] Since in the sequel, we will only be using gamma functions for r < 3, we will write for ri 1 ^, 
F e for ri 2 \ and T+ for ri 3 ', and similarly (with a subscript h) for the hyperbolic versions; we will also omit the 
superscript (r) on P, Q, and R. 

Corollary 3.1. For any parameters fi, v, and any real number x, we have, as v — > + , the estimates (uniform 
over compact subsets avoiding the poles) 

T e (e(vfi + x), e(w - x);e{vuji), e(vu 2 )) _ + v _ ~ CJl ~ ^ m x \\\ 

e(R(v[i; vu>i, vu/2) + R(w; vu/i,vu>2)) 2vio\io2 

-7^7 » r 7 7 ; — r~? w = O^Kt, 

e[R[Vfj,; vu>i, vu>2))i-e{e(yv + x); e{vuii),e(vui2)) 2vlo\U2 

_ °( e ( r^— n<n))), 



T e (e(v/j. + x), e(vv — x); e(vu>i), e(vu} 2 )) 2vuj\u>2 
where $(x) is the continuous, even, periodic function defined by 

§{x) = {x}{l-{x}) = {x}{-x}. 

Similarly, 

e(-R(v(i; vu))0(e(vfi + x);e{vu)) = 0(e(d{x)/2vuj)), 
and likewise for the reciprocal (avoiding the poles). 

Remark. Note that in the above bounds, we can ignore 0(1) terms in R, and may thus replace 

R(vfi;vuj) 1— » — l/12uo; 
R(v[i; vuji , VUJ2) !— > (wi + UJ2 — 2fi)/2AvuJiuj2- 

Also, although we assume x real, the above estimates are clearly still valid if we make a 0(v) perturbation to 
x on the left-hand sides. 

We will thus require some inequalities involving this quantity 
Lemma 3.2. For any sequences c\,. . . ,c n , d\,. . . ,d n of real numbers, we have the inequality 

1<2,j <n 1^*<J^^ l<z<n 
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with equality if and only if the sequences interlace in K./Z; that is, iff they can be permuted so that either 

{ci} < {rfj < {c 2 } < • • < {d„_i} < {c„} < {d n } 

or 

{di} < { Cl } < {d 2 } <■■■< {c n _x} < {d n } < {c„}. 

Proof. First, observe that if two elements Cj, dj agree modulo Z, then their contributions to the inequality 
cancel, and the result thus follows by induction. We may therefore assume that Cj ^ dj for all 1 < i, j < n. 
Now, consider the asymptotics as v — ► + of the case r = w/2 of the determinant identity [7]: 

det e(-l/12Hfl(e(^ + c;^;);eH) 
i<i,i<« ^(e(ur), e(c^ - dj); e(vu)) 

1 > e(n/12voj)6(e(vT);e(vu)) IIi<ij<» %W ~ ^ « 

where real constants (c' ; — Ci)/v, {d[ — di)/v are chosen so that the 2n quantities c^, are all distinct for 
sufficiently small v. Now, since Ci =/= dj, we have 



e(-l/12vu)9(e(vuj/2 + d - d'A:e(vu)) 1 

In particular, the determinant converges to a well-defined limit. Moreover, this limit is nonzero iff the sequences 
interlace, as follows by considering the rescaled determinant deti<j.j<„(sgn({ci} — {dj})). (Indeed, if the 
sequences fail to interlace, two rows or columns must agree; otherwise, the n distinct rows are easily verified to 
be linearly independent.) 

On the other hand, we have the estimates 

e{e(vLj/2 + x);e(vu)) 
9{e(vLu/2);e(vuj)) v v w/ u 

e(l/12vu))6(e(x);e(vu))) = @(e(tf(x)/2voj)) assuming |{x}|, |{1 - x}\ = fi(«), 
giving the estimate 

e(e((#(22ci-di)+ Mci-Cj) + #(di-dj))- ^c l -dj))/2vuo)) 

i l<i<j<n l<zj<n 

for the right-hand side. Since 1/uj) > 0, this is bounded as v — > + iff 

i l<i<_7<n l<i,j<n 

and is bounded away from iff equality holds. The result follows. □ 

Remark. More precise asymptotic calculations give the following well-known (and easy) identity (valid for n > 1) 
as a limit of the elliptic Cauchy determinant: 

det (sgn(x< - Vj )) = 2 "- 1 (_l)«("- 1 ) IT S gn(x. i - Vj ) TT sgn(a; 4 - Xj ) sgn(y 4 - Vj ) 

l<ij<n 



17 



for interlacing sequences with distinct elements, and otherwise. Similarly, for more general values of r, we 
obtain the identity 

det (^n(^)) = (-l)^-i)/2 (x _ 1/xT -i sgn( y c ._ d . )a;Sgn (E < ni<Kj<„ s § n (^ ~ J) ~ d j} 

i<i,j<n if ni<ij<„ sgn(ci - dj) 

We can also obtain a version with hyperoctahedral symmetry. 

Lemma 3.3. For any sequences cq,. . . ,c n , d\,. . . ,d n of real numbers, we have the inequality 

J2 ± dj) - J2 v{c t ± cj) - ± d o) - J2 di ? d ^ ^ °> 

0<i<n,l<j<n 0<i<j<n l<i<j<.n l<i<n 

with equality iff the sequences can be permuted so that 

min({±Co}) < min({±di}) < min({±ci}) < • •• < min({±c n _i}) < min({±<i n }) < min({±c„}). 
Here ± y) := ~d(x + y) + — y) and min({±a;}) := min({.T}, {— x}). 
Proof. Apply the preceding lemma to the sequences ±Ci and 0,±di, 1/2, and use the identity 

§{2x) = 2{d{x) + tf(x + 1/2) - tf(l/2)). 
The given equality condition simply restates the condition that ±Ci and 0,±di, 1/2 interlace. □ 
Corollary 3.4. For any integer n > 1 and any sequence e%, ■ ■ ■ , e n of real numbers, we have the inequality 

2 ^(ei ± ei ) - (n - 1) Y ^(2e. ( )>0 

l<i<j<n l<z<n 

with equality iff the sequence min({±ei}) is constant. 

Proof. The case Co = 0, n = 1 of Lemma 13.31 implies that 

${x ±y) + 2(i?0) - - d{2x) > 

with equality iff min({±:r}) < min({±y}). Adding all specializations of the form [x,y) <— > (ei,ej) with i =/= j 
gives the desired result. □ 

If we rescale Ci,di i— > vci,vdi and take v — > + , the fact that d(vx) = v\x\ — v 2 x 2 for sufficiently small v 
gives us the following limit. 

Lemma 3.5. For any sequences c\, . . . , c n , d\, . . . , d n , of real numbers, we have the following inequality: 
\c%-dj\- Y \ c i-°j\- \ d i- d j\>\ c i~ °M> 

l<z<n l<i<j<n l<i<j<n l<i<n l<i<n 

l<j<n 

with equality iff the sequences can be permuted so that 

c\ < d x < ■ ■ ■ < c n < d n 



d\ < ci < • • • < d n < c n 
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In particular, we have the following fact. 
Lemma 3.6. For any sequences Co, . . . , c„, d\, . . . , d n , of real numbers, we have the following inequality: 

\ci-dj\- Y \ c *- c i\- \ d i- d j\>°i 

0<i<n 0<i<j<n l<i<j<n 

l<j<n 

with equality iff the sequences can be permuted so that 

co < di < ci < • • • < c„_i < d n < c n . 

Proof. Choose a number do such that e?o < mm(co, . . . , c n , d%, . . . , d n , J2i c i ~ J2i>o ^»)> an d a Pply Lemma 
1331 □ 



Remark. The case n = 1 is of course just the usual triangle inequality in K, thus justifying the title of this 
section. 

4 Hyperbolic limits 

Using the above asymptotic estimates for the hyperbolic and elliptic gamma functions, we can obtain corre- 
sponding estimates for the various elliptic hypergeometric integrals of [16] in the hyperbolic limit. In particular, 
in each case, it will turn out that up to an explicit exponential factor, the elliptic integral converges exponentially 
quickly to the hyperbolic integral. 

Let us first consider the case of the Type I (perhaps better named "elliptic Dixon" , see Corollary 17.31 
below and [6]) integral with BC n symmetry, defined for all nonnegative integers m, n, and parameters p, q, 
"o- • • u 2m +2n+3 satisfying 

0<H,M<1, [] u r ={pq) m+1 

0<r<2m+2n+3 

by the integral 

r (m) / x (p;P) n {q;q) n f Yll<i<nYln<r<2 m +2n+?, T e < < U r Z t 1 '-'Pi ( l) T-T dZi 



rii^xn T e {z^zf;p, q) Ui<i< n UK^P, q) ^ n 2tV=1^ ' 



where the contour is chosen to contain all points of the form p l q J u r , < and exclude their reciprocals. 
In the hyperbolic limit p, q, u r — ► 1, this gives rise to the following limit. 

Theorem 4.1. Let /iq, fix,. . . , /J- 2 m+2n+3, ^1, ^2 be parameters such that 

3(cj 2 ) > 0, fi r = (m + l)(cJi + u) 2 ). 

r 

Then as v — > + , 

e(— 2nY^ R(vfj, r ; vto\ , V0J2) + (2n 2 + n)R(0; bui,^))^, ( e ( u A*o) ; e(vfj,i), . . . ; e(vu>i), e(vu>2)) 
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converges uniformly exponentially (over compact subsets) to 

1 f T\l<i<nT\a<r<2m+2n+Z T h{^r ± ^ii^l.^) yr 

(V-wiw 2 )™2™n! J C n ni<i<i<r t r ft( ± ^ ±a; j; w i' w 2)rii<: i <„ r ^( ± 2^;wi,W2) 

where the contour agrees with R outside a compact set, and is chosen to contain all points of the form \i r 
iuji + ju>2, i,j > and exclude their negatives. 

Proof. We first observe that 

e(R(vz; vu>x, VLO2)) 



i(R(0;vuJi,vuj2))(e(vuJi);e(vuJi))(e(vuj2);e(vuj2)) = lim 



z^o (1 - e(uz))r e (e(wz); e(vwi), e(vuj 2 )) 
lim 



(1 - e(vz))T h (z;uJi, w 2 ) 
1 



with uniform exponentially small relative error as v — ► + . 

For the remaining factors, we first assume that Sy(/x r ) > for all r, and thus the elliptic contour may be 
taken to be the unit circle. Now, in the elliptic integral, introduce the change of variables Zj = e(xi), and thus 
dzi/2TTy/—lZi — dxi] this replaces the unit circle by the cube [—1/2, 1/2]". We next claim that if we restrict to 
the smaller cube [—1/4, 1/4]™, the resulting error is uniformly exponentially small. Indeed, we can use Corollary 
13.11 to bound the integrand on the full cube. The \x factors satisfy 

e(-2R(vfJ, r ; voJi,vco 2 ))T e (e(vfj, r ± Xi); e(«wi),e(j)w 2 )) = 0(e( — i- ^ 2 + d{xi))) 

2vLO\L)2 

and thus, using the balancing condition, 

W-(p.(i)ii.~ -\- nr.; V PAiHtU \ p.Hffjoll —Cl(p(- 

2vujiuj 2 



e(-2R(vfi r ; vui,vu2))T e (e(vfir ± £»); e(vwi), e(voj2)) = Q(e( — ■ ^ ^ — 2(n + l)i?(a:j))). 

0<r<2m+2n+3 



Similarly, the remaining univariate factors satisfy 

e(2R{0;vu; 1 ,vu J 2))r e (e{±2x l );e(vLJ 1 ),e(vLJ2)r 1 = 0(e(— (-0(2a*)))) 

2VU>xOJ2 



and, for i < j, the cross factors satisfy 



e(4R(0; uwi, ww2))r e (e(±Xi ± i,);e(wi),e(w 2 )) 1 = 0(e(— - — — (— i?(x* ± a;,-)))) 

2ucjio;2 

Combining these bounds, we find that the integrand is uniformly 

Q(e( ~2?^? ( ^ ( 2 »+ 2 )^)- E ^)- E 0te±*i)))) 

l<i<n l<i<n !<*<.7< n 



Since 



the bound is maximized when 



3( - " 2 ) = 9f(-lM) + 9(-l/wa) > 0, 
ui\uj2 



E (2n + 2)^)- E ^i)- E 
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is minimized, which in turn happens when x\ = x% = ■ ■ ■ = x n = 0, by Lemma 13.31 applied to the case Cj = 0, 
dt = x%. In particular, the integrand is exponentially small everywhere else, and thus restricting to \x%\ < 1/4 
introduces an exponentially small error. 

At this point, using Proposition 12.101 allows us to replace the gamma functions in the integrand with 
hyperbolic gamma functions (times an exponential factor that turns out to be trivial). The factor v~ n from 
((p;p)(q; q)) n can be absorbed in rescaling the variables of integration; we thus obtain the restriction of the 
desired integral to the cube [— l/4w, l/4u]. But again we can bound the integrand, this time using Corollary 
12.31 and find the uniform bound 

Q(e( Z 2^T g ( £ ( 2 ™ + 2)M- E l 2 ^!- E (N + Zil+l^-^l)))), 

1 2 \l<i<™ l<i<n l<i<j<n J 

so the omitted tail is again uniformly exponentially small. 

For the general case, we note that if C is a valid choice of contour for the hyperbolic integral, then for 
sufficiently small v, the image of the subcontour [— l/2w, l/2v] under x i— > e[yx) is a valid choice of contour for 
the elliptic integral. The result agrees with the unit circle outside a neighborhood of size 0(v) of 1; as a result, 
the difference from the unit circle has no effect on the asymptotics. □ 

If we denote the above integral by I^q . h , we have the following corollary, obtained as the limit of the 
corresponding identity for the elliptic case; note that we do not need to compare the exponential factors on 
both sides, since both sides must agree throughout and have generically nonzero limits. 

Corollary 4.2. Let /iq, fii,. . . , ^2m+2n+3 , <^i, ^>2 be parameters such that 

Sj(w 2 ) > 0, y r = (m + l)(^i + u 2 ). 

r 

Then 

l{ BC n ;h(- •• ! Mr,---;Wl,CJ 2 ) = ]J Th ^ r + ^ s ; Wl ' W2 )4c m ;h ■ ■ > "V^ ~ Mr, • ■ ■ 5^1, W 2 ), 

0<r<s<2m+2n+3 

and in particular 

I BC n ;h(--->l J 'r,---;Ui,UJ2) = ]~[ Ifc(/X r + /i s ; Wi, W 2 ). 

0<r<s<2n+3 

Remark. As van Diejen and Spiridonov [5] observed for the Type II evaluation, one can also prove hyperbolic 
results by simply replacing the arguments of [TB] by appropriate limits, rather than taking limits directly. Those 
arguments depend strongly on the fact that the set p z g z generically has finite limit points (in fact, is generically 
dense), which makes analytic continuation trivial. In the hyperbolic setting, the corresponding set Zuji + Zw 2 
is never dense, and only has a limit point when u>i/u>2 is real irrational, so an additional analytic continuation 
argument is needed to extend to generic moduli. 

A similar argument will work in the other cases; some technical issues do arise, however, so it is worth 
discussing those cases as well. 
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For the Type II (again, the name "elliptic Selberg" might be better) integral, the main complication is that 
without an additional condition on the parameters, the integrand is not maximized near Zi = 1. We have the 
following result. Define a family of integrals 

n BC n ( u 0,Ul, ■ ■ • ,U2m+5;t;p, q) 

{p;p) n (q;q) n T e (t;p,q) n f T e {tz?zf;p,g) „ rio<r<2 m +5 T^zf^p, q) dz t 

11 W ±l±l. x 11 



2 "" ! Jo i<£J<„ r e (zf l zf;p, q) x ll n r e (z 4 ±2 ;p, g) W=l* ' 

on the domain t 2n ~ 2 J\ r u r — (pq) m+1 , < \q\, \t\ < 1, where the contour C satisfies C = C -1 , and for all 
h j > 0, contains the points p l q 3 u r as well as the contour p l qHC. 

Theorem 4.3. Let fig, fj,±,. . . , fJ.2m+5, t, u>\, ll>2 be parameters such that 

3(t), 9f(wi), 9f(wa) > 0, (2n - 2)t + ^ ^ r = (m + l)( Wl + w 2 ), 

r 

and satisfying the convergence condition 

-(n-l)r -^-^ 

TTien as t> — > + , 

e(— 2n R(vp r ; vu)\ , VUJ2) + 2n 2 i?(0; vu}\ , uw 2 ) — 2n(n — l)iJ(ur; uwi, W1/2)) 

r 

•#scl( e ( u Mo), e(«/*i), . • . ; e(vr); e(wwi), e(w 2 )) 

converges uniformly exponentially (over compact subsets) to 

I\(r; wi, w 2 )™ /" tt ^(rzfcxjifcxj;^!,^) tt Ilo<r<2m+5 r /t(^r ± x,;a;i,a;2) 

11 rJ+T: + T,-M, r.io'i 11 T\M-9^,-m-, Xi ' 



(^ZJ^) n 2™n!j c „ 11 r ft (±Xi ± xy, uji, lj 2 ) I* T h (±2xi]UJi,u!2) 

l<t<j<n J l<i<n 

where the contour C — —C agrees with R outside a compact set and for all i,j > contains the points 
iu)\ + ju>2 + p r o,s well as the contour iu>i + ju)2 + t + C . 

Proof. Again we change variables to Zi = e(xj) and integrate over the cube [—1/2, 1/2]™; we may also freely 
assume n > 1, as the case n — 1 has already been dealt with. In this case, we find that the integrand is 
uniformly bounded by 

0(e( ~ UJ \~ UJ2 (2 £ #{ Xi ± Xj )-{n-l) J2 0(2**: 

v y 1 A l<i<j<n l<i<n 



-(n — l)r — wi — o^2 



(2(n-l) E ^±^))))- 



(n - l)uwia;2 

By Corollary 13.41 the first $ sum is > 0, with equality iff the sequence is constant; by the case di = of 
Lemma 13.31 the second 1? sum is > 0, with equality iff at most one of the Xi is nonzero. It follows that the 
integrand is exponentially small unless both conditions are satisfied; i.e., unless Xi = 0. The remainder of the 
proof is as above. □ 
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Remark. Note that it also follows from the above proof that the convergence condition is necessary for the 
integrand to be localized. One can readily arrange for equality to hold in the first sum, but not the second, 
at which point if 3(— ((n — l)r + u)\ + u 2 )/wiu;) < 0, the integrand is exponentially larger than its value near 
Xl = 0. 

Let H^BCn-hi^Vi ■ ■ ■ ' T ' w i» W2 ) denote the above hyperbolic integral, as a mcromorphic function on the domain 

3(r), 9f( Wl ), 5(o> 2 ),S( ~ (n ~ 1)T ~ CJl ~ a;2 ) > 0, (2n - 2)r + V Mr = (m + l)( Wl + w 2 ). 

Wiu> 2 — ' 

r 

Corollary 4.4. Lei Mi;- ■ ■ , M5> T , <^i, be parameters such that 

3(t), 3(o) 2 ) > 0, (2n - 2)r + ^/i r = wi +w 2 , 

r 

and satisfying the convergence condition 

UJlLU 2 

Then 

n BC n ;h(^o, ■ ■ ■ ,M5;t;wi,w 2 ) = J| r h ((i + l)r; wi , w 2 ) ]J r /l (ir + ^ r + ^ s ;wi,w 2 ). 

0<i<n 0<r<s<5 

Corollary 4.5. For parameters fio, . . . ,M7; r, Wi, w 2 such that there exists an integer n (necessarily unique) 
with 

(2n + 2)r + ^ fi r = 2(wi + w 2 ) 

r 

and 

9(r),9(^),9(^),9( - (n - 1)T - t,;i -" a ) > 0, 

define 

n h (fj, , ■ ■ ■ , HT, t; wi, w 2 ) := ( ] J T+ (r + /v + r, wi, w 2 ))//j l 1) (r/2 + /x , . . . , r/2 + ^ 7 ; r; Wi, w 2 ). 

0<r<s<7 

TVien 1/^ is invariant under the natural action of the Weyl group E7; in other words, it satisfies the identities 

IlhifJ-o, ■ ■ . ,M7;t;wi,w 2 ) = Ilhivo + v, ■ ■ ■ ,M3 + v,m - v, . . . ,/z 7 - v;t;u\,ui) 
where v = (fi 4 + /i 5 + ^ 6 + ^7 - /x - m - /i 2 - Ai 3 )/4; 

IIh(vo, ■ ■ ■ ,M7;t;wi,o; 2 ) = - /i , ■ ■ ■ ,v- (i, 3 ,i/ - m, . . . ,v' - ht,t;ui,w 2 ), 
where v = (/i + Mi + M2 + M3)/2, = (M4 + /15 + M6 + Mr)/2; anrf 

#fc(Mo, • • •,M7;t;wi,u;2) = /^(zv - Mo, • • • , ^ - Mr; t;wi, w 2 ), 
w/iere f = (/Uo + Mi + ' ' ' + M7)/2; as well as invariance under permutations of /j,q through [17. 
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Remark. Similarly, the other double coset of £7 in Eg that gives rise to (dimension-altering) transformations 
of the elliptic integral also gives rise to transformations of the hyperbolic integral; we omit the obvious details. 
The key observation is that the overall exponential factor that arises when taking the limit is, once one solves 
for n, a function of Yli^h an d is t nus -Eg-invariant. The work of [TS] on recurrences also descends to the 
hyperbolic case; in particular, for t = u>2, one obtains a tau- function for a hyperbolic analogue of the elliptic 
Painleve equation. 

For the A n integral, the difficulty is that the elliptic integral has a condition Yii z i = 1> which in Xi 
coordinates, becomes Xi € Z; this introduces extra complications when maximizing the integrand. Recall 
that the A n integral is defined by 

Ia?{ u 0, ■ ■ ■ Mm+n+l^O, ■ ■ ■ V m+n+1 ]p, q) 

= {p;p) n (q; q) n I IIo<i<n IIo<r<m+»+2 T e (u r z t , V r / Zj ; p, q) dZj 

for < \p\, \q\ < 1, < |u |, . . . , |u m +„+i|, |t> |, ■ ■ ■ , |w m+n+ i| < 1, EL UiVi = (pq) m+1 . (It follows from general 
principles that this can be extended to a meromorphic function on the domain < \p\, \q\ < 1, ]X itjU, = 
(pq) rn+1 , but the condition on the contour is complicated.) 

Theorem 4.6. Let fio, fJ-i,- ■ ■ , fJ-m+n+i, vq, v\,. . . , v m +n+i, u\, ^2 be parameters in the upper half-plane such 
that 

} j fi r + v r = (to + l)(cJi + OJ2) 

r 

Then as v — > + . 

e(— (n + 1) ~S~](R(v fi r ;vcoi,voj2) + R(vv r \ vuji, w 2 )) + {n 2 + 2n)i?(0; uwi,dw 2 )) 

r 

'I A? ( e 0/-*o), e{vfii), e(w ), e(wi), . . . ; e(vw\), e(vuj 2 )) 
converges uniformly exponentially (over compact subsets) to 

Ilo< 4 <« rio<r<m+™+2 F h Or + x i , - Xi J Wl , UJ 2 ) 



f ilo< t <„ ilo<r<m+n+2 Wr + X U ^ ~ Xj^U^ -j-j- ^ 



(V-Wiw 2 )"(n + 1)! J^g.^x^o IL^x™ 1 ^ -^i^j -^;wi,w 2 ) 

Proof. If we perform the change of variables Zj = e(xi) in the elliptic integral, the result is an integral over the 
domain 

-1/2 < x ,Xi,X2, ■ ■ .,x n < 1/2; ^ x i^ z -> 

0<i<n 

a disjoint union of polytopes. Over the entire cube, we find that the integrand is uniformly 

1 2 0<i<n 0<i<j<n 



Now, we find from the case di = of Lemma 13.21 that 



0<i<n 0<i<j<n i 
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with equality iff xq, . . . , x n interlaces with 0, . . . , and Xi € Z; i.e., iff Xi = 0. We thus conclude that the 
integral over the polytope 

-1/A<x ,xi,...,x n < 1/4; J2 Xi=0 

i 

is uniformly exponentially close to the original integral. Thus, as above, the theorem reduces to showing that 
the hyperbolic integral decays exponentially. This in turn reduces to the identity 

(n + 1) \ Xi \- \ x i- x i\ ^° 

0<i<n 0<i<j<n 

with equality only when x± = x% = • • • x n = 0. □ 

The remaining issue in degenerating [16j to the hyperbolic level is the degeneration of the biorthogonal 
functions constructed there. The primary difficulty is that the construction of those functions in [16] does not 
give rise to good uniform asymptotics. However, we can still establish the following. 

Theorem 4.7. Let the parameters To, T\, T2, t%, [1q, fi\, r, u>x, 0J2 be parameters with t, u>i, 0J2 in the upper 
half-plane such that 

(2n - 2)r + tq + t\ + T2 + T3 + (j,q + fix = U\ + cj 2 - 
Then for any partition pair \, and for generic values of the parameters, the biorthogonal function 

U^i- ■ ■ > e ( x i)> ■ ■ ■ i e(vT ):e(vTi), e(«r 2 ), e{vr 3 ); e(w/i ), e(vfj,i); e(ur); e{vuix), e(w 2 )) 

is uniformly bounded for (xi, . . . , x n ) € D(v) n , where D(v) is a region of the form —1 — vC\ < ^R(x) < vC\, 
|^(^)| < vCi, and excluding a hole of radius ev around every pole of the biorthogonal function. Moreover, there 
exists a function 

fc\-hi- --,x t ,.. .;T ;T l ,T 2 ,T 3 ;iM ) ,fj, 1 ;T;cJi,U2) 

such that as v — ► + ; 

Hx\- ■ ■ , e(xi), . . . ; e(vT ):e(vTi), e(vn), e(ur 3 ); e(v/J,o), e(vfj,i); e(vr)\ e(vwi), e(vu2)) 
~ K\ih(- ■ -i x i/ v > ■ ■ ■ \ to:ti,T2,t 3 ; (io, Hi;t;u}i,W2) 

converges exponentially to 0, uniformly for x in a compact subset of the domain — 1 < 5R(a;) < 1 . 

Proof. We first observe that the claims of the theorem are certainly true if we replace TZ^ 1 by a product of 
functions of the form 

or similarly for k>2- I n particular, it was established in [16) that there exist functions iv" (jjlo:;t',Wx,W2) of the 
above form such that there exists an expansion 

^ = E cw* 

/iCA 
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for some coefficients C\fj, independent of z$. It thus remains only to show that for generic parameters, these 
coefficients Ca m converge exponentially. Moreover, the action of the integral operators of [16] can be computed 
explicitly in the F\ basis, and the coefficients of the corresponding matrices converge exponentially (to a 
triangular matrix with generically nonzero diagonal). Thus the generalized eigenvalue equations satisfied by 

~ (n) 

K^> set up linear equations in the Ca m with exponentially converging coefficients. Since the limits of the 
generalized eigenvalues are generically distinct, the limiting linear equations are generically nonsingular, and 
the result follows. □ 

Remark. The unviariate hyperbolic biorthogonal function 

^A;A ( X ! T :7 "1 ) T 2 J T 3 ! MO i Ml I T\ U>1 , L0 2 ) 

was discussed in §8.3]. 

Note in particular that if we multiply the integrand of either BC n integral by a function satisfying such 
convergence properties, the resulting integral will also converge exponentially (assuming the unadorned integral 
so converges). Also, a similar argument works for the interpolation functions (which as special cases of the 
biorthogonal functions do not quite fall under the above generic result, but again satisfy suitable integral 
equations). As a result, every identity of [16) involving such functions converges exponentially (possibly with an 
explicit factor of the form exjp(a + bv)) to a corresponding hyperbolic limit. One should note (as observed in [201 
§8.3]) that further degeneration of the parameters can lead to convergence issues, as without the moderating 
effect of the poles, the biorthogonal functions grow exponentially as |3?(x)| — > oo. 



5 Trigonometric limits 

The main difficulty with the trigonometric limit p — > is that the general case of the transformations involves 
parameters tending to infinity, making the contour ill-behaved in the limit. This can be fixed at the expense of 
breaking the symmetry of the integrand. 

Recall that for the type I BC n integral, the parameters are constrained to satisfy the balancing condition 

IJ u r = (pq) m+1 . 

0<r<2n+2m+3 

The natural way to satisfy this in the p — > limit is for 2n + m + 3 of the parameters to be 0(1), while the 
remaining m + 1 parameters are Q(p). This then makes the p — > limit of the integral trivial to compute. 

Theorem 5.1. For any parameters uq, . . . , U2 n +m+2> u 0> • • • > v m, 9 satisfying 

\q\ < 1, IJ U r = Y[ V r, 

0<r<2n+m+2 0<r<m 

we have the limit 

Ji^sct^o, ■ ■ .,u 2n +m+2,pq/vo, . . . ,pqjv m ;p, q) 



2 n n\ 



(<7; ( ?)" f TT w ±1 ±1 s-l TT no<r<2n+m+2 r t( M '' Z l ±1 ;9) dZi 
J cn 11 Z J 5 9) 11 w „±2.^tt FT! .±1. 



l<i<j <n 



ii<n V ^ Z ? ; i) Uo< r < m rt(v r z?\ q) 2n^l, 



2G 



where the contour contains all points of the form p l q : >u r , i,j > 0, and excludes their reciprocals. 
Proof. This follows immediately from the facts that as p — > 0, 

r e (a;;p,g) ±1 = T t (x; q)^ (1 + 0{p)) 
T^pq/x^q)^ = T t (x;q)* 1 (l + 0(p)), 

uniformly in x away from the poles. □ 

Unfortunately, the right-hand side of the type I transformation involves parameters 

(pq) 1/2 /u , (pq) 1/2 /u 2n +m+2, (pq) ~ 1/2 v , (pqy 1/2 v m , 

which as mentioned above gives an apparently ill-behaved limit. The primary difficulty is that the divergent 
parameters not only deform the contour, but in fact pinch the contour in the limit, making it approach both 
and infinity. It turns out, however, that there is a way to break the symmetry in such a way as to eliminate 
half of the offending poles, thus allowing the contour to be renormalized, giving a well-behaved limit. 

The key fact is the following identity of (/-elliptic functions. Here R{z{) denotes the operator such that 
R(zi)f( Zi ) = f(l/zi). 

Lemma 5.2. For any parameters u , . . . , u n+ \, q we have the identity 

TT , o, . gtlWl Ur/ IIl<,<» fij g) IWn gKfjjg) TT ft( , 

11 + — ff e(z'z~a) " " = 11 6{u r u s -q). 

<s<n+l 

Proof. The left-hand side can be expressed as a sum of 2™ terms, all of which are elliptic functions in z with 
respect to multiplication by q, and thus the sum is also an elliptic function. Moreover, since the original function 
is invariant under permutations, the sum is invariant under the action of BC n . In particular, the order of the 
sum along each reflection hypcrplane must be even; since the summands have at most simple poles there, it 
follows that the sum is constant. The constant can be recovered by taking zi = Ui, making all but one summand 
vanish. □ 

Lemma 5.3. For any nonzero parameters to, ■ ■ ■ , t n , Uo, ■ ■ ■ , u n + m +i, vo, ■ ■ ■ j v m, P, q with 

o<M,M<i, Y[ tr= EI Ur 1 1 Vr 

0<r<n 0<r<n+m+l 0<r<m 

and any complex parameters a,w/0, we have the identity 

Yl 9(t r t s I 'a; q)I { ^ n (t / a 1/2 , . . . , t„/a 1/2 , a 1/2 /u , ■ ■ ■ , a 1/2 /u m+n+1 ,pq/a 1/2 v , . . . 1 pq/a 1/2 v m ;p,q) 

0<r<s<n 

{p;p) n {q;q) n f ^(rio<r<„ t '-/wni< J <„^;g)ni< l < Il ^(^/w;9) Hl<i<j<n 6 (P Z i Z j/ a '>P) 



™ ! Jen Uo<r<n ( t r/^-<l) Ul^j <n T e {{Zi/Zj ) ±X \ p, q) 

nrr r / . / , , x tt r e (zi/u r ;p, g) 

Ki<«0<r< n 0<r<m+„+l ^ e(PqZ t U r / a, p, q ) 



Y e {pqzi/av r ;p,q) dz % 
; r e (z iVr ;p,q) 2ir^T Zl 



n 
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where the contour contains all points of the form p l qH r , p l q 3 a/u r) p % (j? /v r , i,j > 0, and excludes all points 
of the form p~ 1 ~ I q~ J a/t r , p~ l q~3u r , p - " 1 q~ 3 av r , i,j > 0. 

Proof. If we multiply the integrand on the left-hand side by the case u r = a^ x ^ 2 v r , < r < n, u n +i = a 1 ! 2 jw 
of the lemma, the symmetry of the integrand implies that E[i<i<n(^ + R( z i)) can t» e replaced by 2 n . Shifting 
the variables of integration by Zi —> a~ x / 2 Zi gives the right-hand side, up to a shift in contour with no effect on 
the integral. □ 

Remark 1. Note that for specific choices of w, the contour condition may conceivably be weakened; the point 
is that the to-dependent factors can cancel out poles of the integrand, making the corresponding constraints on 
the contour superfluous. In particular, for certain specializations of the parameters, it can be the case that the 
contour conditions for generic w are inconsistent, but a suitable choice of w makes the integral well-defined. 

Remark 2. If we multiply the integrand on the left by a symmetric function / (adjusting the contour conditions 
accordingly), the effect is to multiply the nonsymmetric integrand by /(. . . a~ x / 2 Zi . . . ), with suitable contour 
conditions. 

This makes the limit p — > trivial again, as long as \pq\ < \a\ < 1. Taking a = pq gives the following. 
Theorem 5.4. For any nonzero parameters uq, . . . , U2n+m+2, vq, . . . , v m , q with 

0< |<7| < 1, J] U r = Yl V r, 

0<r<2n+rn+2 0<r<m 

we have the limit 

lim n Fl d ( v r v s/pq;q)lBc((pq) 1/2 /u 0l • • • , (pq) 1/2 /u 2n + m +2 1 (pq)~ 1/2 v , . . . , (pq)~ 1/2 v m ;p,q) 

Q<r<s<m 

= (g; q) m r ^(IIo<r<m ^vMTi^m z ^ g) IIi<i<m 0{zi/w, g) rii<i< 3 < m ( 1 - gjj/gj 

m! Jc» Uo<r<m e ( V r/ w >l) Ul<i<j<m T t{{Zi/ Z,)^] q) 

n n Ur, \ n r * k^/ g, v r / dz ^= - , 

l<i<m0<r<2n+m+2 l t\ Z i U r^) Q < r < m M^J-lZi 

where the contour contains all points of the form q 3 V r> j > 0, and excludes all points of the form q~ 3 u r , q 3 /v r) 
j>0. 

Corollary 5.5. The trigonometric integral of Theorem \5.1\ is equal to 

Y[ T t (u r u s ;q) J| T t (v s /u r ; g) -1 J][ T t (v r v s /q; 

0<r<s<2n+m+2 0<r<2n+m+2 0<r<s<m 

0<s<m 

times the trigonometric integral of Theorem \5.4\ 

Remark. The univariate cases n = 1, m = and n = 0, m = 1 are the Nasrallah- Rahman integral and an 
integral identity of Gasper (equations (6.4.1) and (4.11.4) of [5]); the general m = case is due to Gustafson 
i- 
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We also obtain a nontrivial transformation by taking a ~ p a for < a < 1, say a = (pq) 1 ^ 2 for symmetry. 
Theorem 5.6. For any nonzero parameters to, ... , t n , uq, . . . , u n+m+ i, vq, ■ ■ ■ , v m , q with 

o < \q\ < i, n ^ = n w r j| u r 

0<r<n 0<r<n+m+l 0<r<m 

and any complex parameter w, we have the limit 

hm II 0{(pq) 1 t r t s ;q)r BC !{ ... , ;p,q) 

P 0<r< S <n ^ U ° u m+n+l V V m 

{q;q) n f ^I\ G <r<n t rlwY\ 1<l<n z l ;q)Y\ 1<t<n e{z l /w-q) „ 
n\ Jc~ n <r<nWr/w,q) 

TT IIo<r<K ^tjtr/Zj] q) IIo<r<m+n+l ^t(zj/u r ] q) dZj 
^ Uo<r< m r t( Z i V r;<l) 2-K^lz, 



l<i<n 

where the contour contains all points of the form q l t r , i > 0, and excludes all points of the form q~ l u r , i > 0. 
Corollary 5.7. The trigonometric integral of Theorem \5.6\ is independent of w, and if multiplied by 

n n utr/u s - q )-\ 

0<r<n 0<s<n+m+l 

is invariant under the involution 



{ji, m, . . . j v r u r t r , . . . ). 



Remark. This can also be obtained as a limit of Corollary 15.51 after first breaking the symmetry of the left- 
hand side as in Lemma 15.31 In particular, this should perhaps be thought of as a degeneration rather than a 
direct limit; we mention it to point out that that distinction is somewhat artificial (any degeneration should be 
obtainable as a limit directly from the elliptic level), but more importantly because the Type II analogue has 
important consequences. 

For the Type II integral, we again have a trivial limit in one case. 
Theorem 5.8. For any parameters uq, . . . , u TO +4, Vo, ■ ■ • , v m , q satisfying 

\q\<l, t 2 "- 2 [] U r= II «n 

0<r<m+i 0<r<m 

we have the limit 



lim , n BC<( u o> ■ ■ -,Um+i,P<l/vo, ■ ■ ■ ,pq/v m ;t;p, q) 

p—*0 n 



(q;q) n T t (t;q) n f yr ^(tzf 1 zf-q) „ rio<r<m+4 ^(ttr^f 1 ; q) dz t 
11 r.r„±i„ ±i.^ 11 



2 n n\ 



° n l<i<j<n 



Uz^zf-q) ^ <n T t (z± 2 ; q) Uo<r< m U^zf^q) 2ttV=1* ' 



where the contour C satisfies C = C , and for all i > 0, contains the points q l u r as well as the contour q % tC ' . 



29 



We recall from [16] the following identity, which plays the role of Lemma T5. 2 1 for the type II BC n integral. 

n,, , R , ss 9(u a z u Ul Zj, u 2 z ly t^upu^/zi; q) -p-r Q(tz lZ ^q) -p-r 
(l + R{Zi)) e(z 2 -a) H 9(zz -a) = H ( U °" 2 ' Ul " 2 ' 9 ^ 

l<i<n \ i>1) l<i<j<n ^ 1 0<i<n 

Lemma 5.9. For any nonzero parameters tp, t\, uq, . . . , u m +i, vq, . . . , v m , p, q with 

o<ipi,i?i<i, ^ 2 wi= n u r n v r 

0<r<m+2 0<r<m 

and any complex parameters a,!c/0, we have the identity 
II 0(t%ti/a; l) n BC n (to/a 1/2 , ii/a 1/2 , a 1/2 /u , . . . , a 1/2 /w„ i+2 , M /a 1/2 v , . . . ,pq/a^ 2 v m ; t; p, q) 

0<i<n 

_ {p;p) n (q;q) n ^e(t;p 7 q) n f rr 6{zi/w,t n ~Hoti/wZi;q) -p-r T e (tpZjZj / a,pqZiZj / a,t{zi/ Zj)*- 1 -^, q) 
n\ J C n eip-Ho/Wjt^h/wiq) 1< ^ <n T e (pziZj/a,pqziZj/ta,{zi/zj) ±l ;p,q) 

Yl 9(pzf/a;p)r e (pt z i /a,pt 1 z l /a,t /z i ,t 1 /z l ;p,q) 

l<i<n 

T e (zi/u r ;p,q) -p-r T e (pqzi/av r ;p,q) dz l 



n L e yz l /u r ;p, q; -p-r 



o<r<m+2 re(pg^u r /a;p, g) ^.1 to r e (z 2 zv;p, 9 ) 27rV=T* 
where the contour C is chosen so that for all i,j > 0, it contains the points and contours 

pV to, pV'ti , ?Vfl/«r, P i+ V +1 K, p¥tc, pV'ta/c 

and excludes the points and contours 

a/p i+1 qH , a/p i+1 qH 1 ,u r /tfq\av r /tf +1 qJ +1 , a/p i+1 qHC, C/p l qH. 

Remark. It is possible to choose a contour of the given form satisfying C = aC^ 1 , namely a~ 1 Co where Co is 
a suitable contour for the left-hand side. 

Theorem 5.10. For any nonzero parameters uq, . . . ,M2m+3, Vo> vi, q with 

< |g| < 1, Y[ u r = t 2n - 2 v oVl 

0<r<2m+3 

and any complex parameters w ^ 0, we have the identity 

lim TT e^vovilm^Bc (vo/(pq) 1/2 ,v 1 /(pq) 1 / 2 ,(pq) 1 / 2 /u ,...,(pq) 1 / 2 /u 2m+3 ;t; P ,q) 

p— >0 A A ™ 



0<i<n 



(q; q) n T t (t; q) n f -p-r Q(zj/w,t n 1 v v 1 /wz l ; q) y\ ^t{tZiZj /q, ZiZj,t(zi/ 'zj) ±l \q) 



' c "i<i<„^ ^o/w,^ l vi/w;q) i ^ j<kn Tt{z i z ] /q,z l z j /t 1 {z i /z j ) ±1 -q) 



T t (zi/u r ;q) dz, t 



II (1 - 2 2 /g)r 4 (w ^/<?,'fiz l /q,wo/z l ,'i;i/z,;g) J| 

l< t <n 0<r<2m+3 i *^ u r, 9J -JttV-I*. 

where the contour C is chosen so that for all i > 0, it contains the points and contours 

q l v ,q l V!, qHC, 
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and excludes the points and contours 

q^/vo, q^/vi.Ur/q*, q 1 ^ /tC, q^C/t, 

assuming such a contour exists. 

Remark. It is easy to verify that there exist choices of the parameters for which a circular contour of radius 
q 1 / 2 satisfies the given conditions, and thus the integral on the right has a well-defined meromorphic extension 
to general parameters (and the limit will continue to hold); the only question is whether this can be obtained 
from a domain of integration of the form C n . 

Corollary 5.11. When m = 0, the above trigonometric integral evaluates to 

T t (t i+1 ,t l v v 1 /q; q) Ilo<r<3 ^t(t l v /u r , t l v x ju r ; q) 



n 



We also obtain a transformation. 
Corollary 5.12. When m — 1, the trigonometric integral of Theorem \5.8\ is equal to 



n 



U.0<r<s<5 T ( tiu rUs;P, (?) 



,, „ n o <r<5 r ( w o/^Ur, , yi/^M r ;p,g)r(uowi/^g;p,9) 

times the image of the trigonometric integral of Theorem \5.10\ under the specialization Ui ► t^ n ^^ 2 Ui, Vi — ► 

t -(n-l)/2 v ._ 

There are other transformations relating these integrals, but all can be obtained by applying the above 
transformation to one or both sides of a transformation of the integral of Theorem 15.81 alone. 
Similarly, taking a — p 1 / 2 above, we obtain the limit 

Theorem 5.13. For any nonzero parameters to, t\, uq, . . . , u m +\, vq, . . . , v rn , p, q with 

o<h,m<i, n Ur n Vr 

0<r<m+2 0<r<m 

and any complex parameter mi/0, we have the identity 

lira TT ^ftotx/p 1 / 2 ;^//^ (t Q / P 1/4 , Wp V V /4 /«o, • ■ ■ ,P 1/4 /u m+2 ,p 3 ^q/v , . . . ,p 3 / 4 q/v m ;t;p,q) 

0<i<n 

_{q-q) n T t {t-q) n f ^ 0( Zi /w, t n ~ x toh/wz,; q) T t (t( Zi / z^ 1 ; q) 



ftl+i-^+ n /on /i/r 11 

<i< 

n 



r t (* / Zi,ti/Zi] q) ]lo<r<m+2 T t {Zi/u r \ p, q) dz t 



where the contour C is chosen so that for all j > 0, it contains the points and contours 

qftorfh, qHC 

and excludes the points and contours 

Ur/qi, C/qH. 
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Of course, with the Type II integral, we are particularly interested in the effect of multiplying the integrand 
for to = by the biorthogonal functions. Note that since the integral is taken over a compact curve in each 
case, the limiting relation will continue to hold as long as the limit of biorthogonal functions exists, and (more 
difficult) the revised contour conditions are satisfiable in the limit. The primary constraint is that we may 
only consider p-abelian biorthogonal functions, since otherwise the contour must contain at least one point 
converging to oo as p — > 0. For the first two limits, there is no difficulty with convergence of the biorthogonal 
function. Indeed, the p-abelian biorthogonal functions satisfy the further identities 

= ^ ) (.../^...;^:p-A,A ) A;P 1 \p- 1/ Vii; M ) 
= ^o™ (• -- z i-- ■ih-h,h,h;uo,ui;t;p,q) 

and thus in each case the relevant limit of biorthogonal functions is the same. 

For the a = y/p limit, the situation is more delicate, but we find that if t 2n ~ 2 aoaibobicd — q, then we have 
a well-defined limit 

R £as-i(- ■ ■ ' z *> ■ ■ ■ > a o-a-i,b , h; c, d; q, t;p) 

:= li m ^H-.-p- 1/4 %.--;p- 1/ V:p- 1/ Vy /4 6oy /4 6i;p 1/4 c,p 3 / 4 d;«;p,g) 

= lim ..p-^/zi . . . ;p 1 / 4 a : P 1 / 4 a 1 ,p- 1 /4 6o , J3 -i/4 6i;p 3/4 Cjp i/4 d . ^ q) _ 

p— »o 

This is a multivariate analogue of the biorthogonal rational functions of Al-Salam and Ismail pQ. More precisely, 
by specializing the a = ^Jp limit appropriately (in particular, w — c _1 ), we find that the functions 

R x?AS-i(- ■ ■ ' z *' • ■ ■ ;a :ai,6 ,6i;c, d;q,t;p) and R ( ™ AS -i(- ■ -A/zi, ■ ■ ■ ; &o:&i, do, ay, d, c; q, t;p) 



are biorthogonal with respect to the density 



n r t ft(z i /z J -) ±1 ;g) y\ Vt{ao/z il a 1 /z ll b z il b 1 z i ;q) 



Ttiizi/zj)^ 1 ] q) }i T t (q/czi,qzi/d,t n 1 a aic/z i ,t n 1 bobidZi;q)' 

l<i<J<n l<i<n x ' 1 ' 

which becomes Al-Salam and Ismail's density when t n ~ 1 aoaic = q 1 ^ 2 = t n ~ 1 bobid and n = 1. The constraints 
on the contour are independent of A, fx, c and d, and are simply that C must contain the points q l a r and 
the contours q l tC, and exclude the points l/q % b r and the contours C/q l t. If we then take a\,b\ — > and set 
a = b = q 1 / 2 , we obtain polynomials biorthogonal with respect to the density 



n 



,_, ( B r t ((*/*,)±i;,)' 

these are, of course, simply the ordinary Macdonald polynomials up to a suitable normalization. That 
these arise as limits of the biorthogonal functions is not particularly new (since they are limits of Koornwinder 
polynomials); what is new is that a limit exists that respects the inner product. 

It should be possible to obtain similar limits in the A n case; since the contour conditions are significantly 
more complicated in that case, however, we mention only the identity which presumably plays the role of Lemma 
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15.21 in this case: 

U.O< t <n °( xz */ U.O<j<i V 3 Ui<j<n U f> «) U.O<j<i °i v J z ^ O) Ui<j<n d i u ] z u O) 



symm s - 



n <i<j< n ZjO^/z^qjUjO^/uj-q) 

n 2 j ;(? ) n 6 ( x / n ^ n 



= V{X 

0<^<n l<i<n 0<j<i 



a special case of Theorem 4.4 of [PT] . 



6 Rational limits 

The rational limit is most naturally viewed as a combination of the hyperbolic and trigonometric limits, and 
thus in particular requires both the asymptotic calculations from the hyperbolic case and the symmetry breaking 
from the trigonometric case. In addition, it can be reached by taking lu2 — > in the hyperbolic case, or q — ► 1 
in the trigonometric or elliptic cases. We consider the limit from the elliptic level, as the other levels introduce 
no further complications. In each case, the integrand factors as a product of g-theta functions and functions to 
which Corollary 12 . 151 applies: the exponential behaviour of the integrand comes only from the former. 

Theorem 6.1. For 9f(w) > ; Eo<r<2n+m+2 AV - Eo<K m y r, 

( \ n(2n+3) 
(p;p)y2vrucj/V T Tj 

^BC n ( e ( v Mo), • ■ ■ ,e{v^2n+m+2),pe{v{uJ - vq)), ■ ■ ■ ,pe{v{uJ - v m ))\p,e{vu)) 

(V2^Lj/V^l)- n f jj X-l TT n <,.<2»+m+2 r r(^ ±Xj-,Uj) 

= ^Tl / r r (±Xi±Xj;u) II - - - — — rdxt, 

where C is a contour agreeing with M outside a compact set, and separating the points of the form fi r + ju>, 
j > from the points of the form —fi r — ju, j > 0. 

Proof. We consider the case in which the original contour is the unit circle; deformed cases are analogous. As 
in the hyperbolic case, we make the change of variables Zi — e(xi) and integrate over [—1/2, l/2] n . If we divide 
the integrand by 

\{ 9{e{ Xl +x 3 )-e{vuo)) JJ 0(e( Xi - ay); e(vu)) J[ 9(e(vuj/2 + Xi ); e{vLo)y 2n - 2 

l<2<j<n l<i<j<n l<z<n 

the remaining factors of the integrand are controlled by Corollary 12.151 to have at worst polynomial growth in 
v, and thus the exponential behavior of the theta functions dominates. In particular, up to polynomial factors, 
the integrand satisfies the bound 

0(e((-l/2«w)[ 2(n+l)0(x t )- J2 ^ l ±x ] )- ^ 0(2*,)])) 

l<i<n l<i<j<n l<z<n 

which decays exponentially unless x%,... ,x n = o(l). We can thus restrict the integral to [—1/4, 1/4]™ and rescale 
the variables by v. The result then follows from Corollarv l2.14l □ 
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The other case is more complicated, in that the exponential contribution to the asymptotics is not enough 
to properly localize the integral. 

Theorem 6.2. For 9f(w) > 0, Eo<r<2n+m+2 Mr = E <r<m v n 

/ , \ 2m 2 +3m 

lim e(-(m 2 + 2to)/24vcj) I (p\p)\J 2kvuj / ][J 0(e(v(v r + v s - uj))/p; e(vu>)) 

» ' 0<r<s<m 

r(«) (rill-, \\ rill-, \\ e{v{v -u/2)) e(v(v m - u>/2)) 

Vp Vp 
m! 7 cm rii<i< ro ^(^ - «>; w ) _1 n <r<m - w ; w ) Ui< i<3 < m r r(±(^ - ^) 

I! IT ^~T^ , ^ r '^| TT r r (i/ r + af» - oj, v r - Xi;w)dxi, 

Ki<m0<r<2n+m+2 v ^ " / 0<r<m 

where C is a contour agreeing with R outside a compact set, and separating the points of the form v r + jw, 
j > /rom </ie points of the form fi r — juj, —v T — (j — 1)lo, j > 0. 

Proof. We begin with the integral of Lemma 15.31 replacing the extra parameter by e(vw). The exponential 
factor in the asymptotics of the resulting elliptic integrand is 

l<z<m l<2<j<m l<i<m 

and thus the integrand is exponentially small unless the sequence xi,...,x m in R/Z interlaces (or nearly 
interlaces) with the all-zero sequence. More precisely, if we split the integral into 2 m integrals based on the 
decomposition R/Z = [—1/4, 1/4] U [1/4,3/4], then any piece with more than one [1/4,3/4] factor contributes 
an exponentially small amount. For the pieces with exactly one [1/4,3/4] factor, we find that upon rescaling 
the [—1/4,1/4] variables, the resulting integrand has order 0(v 2 ) and is integrable; thus those pieces again 
contribute a negligible amount to the limit. We may thus restrict our attention to [—1/4, 1/4]", or cquivalcntly 
(up to 0(v)), the integral over [— l/4u, l/4u] of the rational limit integrand. The omitted tails are again either 
exponentially small or have integral of order 0(v 2 ), so the result follows. □ 



Corollary 6.3. The rational integral of Theorem \6.1\ is equal to 

Yl T r (n r +n s ;uj) Y[ T r (v s ~ fi^wy 1 J| T r {v r + v s - uj; u>)~ 

0<r<s<2n+m+2 0<r<2n+m+2 0<r<s<m 

< s < m 



times the rational integral of Theorem \6. e A 

For the Type II integral, again the first case is straightforward. 
Theorem 6.4. For any parameters fi , . . . , fi m +4, z/ , . . . , v m> u>, r satisfying 

9f(w),9f(T)>0, (2n-2)r + J] Mr . = J] u r> 

0<r<m+A 0<r<m 
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we have the limit 

2n(2n~3)T/oj+6n 



lim + e(— n/Avuj) ( (p;p) y 2-kvuj/^^I 

T {m) 



#Bc„( e ( w Mo), • • • ,e(v/j, m+4 ),pe(v(u - v Q )), . . . ,pe(v(u - v m ))\ e(vr); p, e{vw)) 

r r (r;w)" f yr r r (r ± X z ± XjjLj) yr I\a<r<m+i^ ril^r ± Xj;U}) 



w/iere i/ie contour C agrees with K outside a compact set, satisfies C = —C and for all i > 0, contains the 
points iuj + \i r as well as the contour iu; + t + C. 

Proof. The exponential factor in the asymptotics of the integrand is 

0(e(%2(#(2xi) - 40( Xi ))/2vu)), 

i 

which is exponentially small unless Xi = 0. The limit follows as above. □ 
The nonsymmctric Type II integral has even worse behavior than the nonsymmetric Type I case, however. 
Theorem 6.5. Let /j, , . . . ,^ 2m +3, v<a, v\, oj, t, p be parameters such that \p\ < 1, 3(oj),3(t) > 0, and 

(2n - 2)r + v + v x = ^ fj, r , 

r 

as well as the convergence condition 3J(t/w) > — 1/n. Then we have the limit 

(I \ 2n(2n— 3)r/w+6n 
{p;p)\J2irvoj/V^l\ ]] 8(e(v(iT + v + v\ — uj)/p); q) 

' Q<i<n 

<( C( ° ( y /2) , e -^-=^ ,^pe{v{u/2 - mo)), • • • , y/pe(v{u>/2 - ^ +3 ));e( W );p, e(vw)) 
T r (T;uj) n / -j-j T r (T + Xi + Xj — u>,Xi + xj,t + Xi — Xj,T + Xj — XijU)) 



(V2iruj/y/^l) n nl Jc" x<i<j<n ^r{xi + Xj -u,x t + xj -T,x t - Xj,Xj -x^uj) 

I [ (2x 1 /lu - l)T r (v + Xi - oj, v x + Xi - uj, v - x i} v x - x^w) 

l<i<n 

9 h (xj - w, (n - 1)t + vp + vx - w - x^lj) yr T r (xj - fj, r ;u) 
0h{{i - 1)t + u -w,(i- 1)t + vx - w;uj) Q<r ^ m+3 ^r{xi + (i r ;u) 

where the contour C is chosen so that for all i > 0, it contains the points and contours 

iu) + vo,iu) + vx, iuj + r + C, 

and excludes the points and contours 

(1 — i)u — i/o, (1 — i)u) — vx, jj, r — iu>, (1 — i)u> — t — C,C — ito — r, 
assuming such a contour exists. 
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Proof. The exponential factors in the asymptotics of the elliptic integrand actually cancel completely, with the 
result that the integrand has polynomial asymptotics. The result will follow from dominated convergence if we 
can show that the rational integrand converges. 

For the rational tails, we find that (assuming C = uj/2 + R for simplicity) the integrand converges iff the 
integral 



converges for e > 0; equivalently, via the change of variables yi = x 2 /e 2 , the theorem reduces to the convergence 



Up to linear fractional transformation, this is an instance of the Selberg integral ([25], stated as Corollary 17.51 



7 Classical limits 

The final limit we consider is that corresponding to the usual beta integral. Although the beta integral itself is 
generally viewed as the bottom level, this is in fact a somewhat misleading view, as the integrals we obtain are 
in fact still elliptic (involving powers of theta functions). For the BC n cases, a suitable change of variables exists 
that essentially eliminates the dependence on p, but the corresponding change of variables for the A n integral is 
much less obvious (if it exists at all). Furthermore, even for the beta case, the classical transformation analogue 
can only easily be reached by degenerating either the hyperbolic or elliptic levels; the symmetry breaking of the 
trigonometric and rational cases introduces unnecessary complications. 

Since the integrand involves powers of theta functions, there are in general some subtle issues involving 
choices of branch. It will thus be convenient to restrict to the case p real, where the phases are easier to control. 
We have the following. 

Lemma 7.1. Choose p and z such that —1 < p < 1 and \z\ = 1. Then the standard branch of log9(z;p) 
satisfies 




of 




below), and thus converges as long as 31(t/u>) > —1/n, as required. 



□ 



log0(z;p) = log(-z)/2 + log \9{z;p)\, 
log6(p 1 / 2 z;p) = log\6(p 1 / 2 z;p)\. 




taking the principal branch of the logarithm on the right-hand side. Now, 



log(l - p l z) + log(l - p l /z) = 2 log |1 - p*z\, 



so it suffices to show that 



log(l -z) = log(-z)/2 + log |1 - z\, 
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which follows from the observation (1 — z)/ 'y 1 —z = \1 — z\. Similarly, 

\og6{p 1 / 2 z-p) = ^21og|l-^ +1/2 z|. 

0<i 



□ 



We will thus assume —1 < p < 1 in the sequel; note, however, that the case of more general p can be 
obtained by replacing 

\9{z;p)\" ^ {-zY' 2 6{z-pY 
\6i P 1/2 z; P )\ K ^6(p 1 ^z;p) K . 

In any event, we need all parameters to have absolute value 1, IpI 1 / 2 , or \p\ within 1 + o(l) for this to work. 

For the Type I BC n integral, it is particularly natural to take 2m + 2 parameters to have norm jpj 1 / 2 , at 
which point both sides of the transformation take the same form. More general cases could be considered, but 
appear to give rise to the same limiting identities, so we will restrict our attention to the simplest case. 

Given points x, y, z on the unit circle with x, z distinct, y € [x, z] or equivalently x < y < z indicates that y 
is on the closed counterclockwise arc from x to z, and similarly for open arcs. 

Theorem 7.2. Let clq, . . . , a n , bo, . . . , b m be points on the unit circle with 

1 < ao < a\ < ■ ■ ■ < a n < —1, 

let to be a point in the upper half-plane, and let a^, . . . , , f3^ , . . . , (3^ be parameters such that 5ft(a^) > and 

< r < n < r < rn 

Then, writing q = e(vtu), a r = a+ -\- a~ y f3 r — /3+ + f3~ } we have 

r r e (q^^;p,q) 

lim T T T = T = 

v ^ 0+ no<r<« Fe(q ar ;p, q) n <r<,<n r e (9 Q - +Qs a r a s , q a r +^ a r /a s , q a ^ + a ° a s /a r , q^ /a r /a s ;p 7 q) 
. .,q<*ta r ,q<*;/a r , ....... /b r ,p^q^-^ b r , ...;p,q) 

0<r< S <n llo<r<n L \ a r ) 

„±l._\|a r -l 



J ZiZ[ai-L<H] iKi^Kn l<i<„ I lO<r<m \°U> Mj ,P)\^ lZ t 

Proof. Using Lemma 13.31 we find that the integral decays exponentially outside the stated product of arcs (or 
images under the hyperoctahedral group). The result then follows by dominated convergence. □ 

If we define 4>{ z ) — S{z;p) 2 /9(—z;p) 2 , then 

9(-l;p) 2 z9(wz ±1 ;p) 



(z) - cj)(w) 
<U{z) 



9(-z;p) 2 9(-w;p) 2 
9(-l;p) 2 (p;p) 2 9(z 2 ;p) 
6{-z;pY 
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so for z in the arc [1,-1], 



\0(-z;p)\ 4 2ttV^Tz 
and 

Consequently, we can make a change of variables in the resulting transformation to obtain the following result 
of Dixon [6] . 

Corollary 7.3. For any parameters ao, . . . , a n , bo, . . . , b m , ao, ■ ■ ■ , «n, Po, ■ ■ ■ , Pm such that 

-b m < -b m -i < ■ ■ ■ < -b < a < oi < • • • < On, < 3?(a r ),3?(/9 r ), 2J a '- = X! ^ r 

0<r<n 0<r<m 



we have the identity 

no< I<3 <ni^-%r~ a '^ r ii < i < n \aj-xi\ a '- 1 
no<i<nrK)n <,-< m i«i+^h i 4 e K- 11 a i ] 1 <y<™ 1 Ji no<i< m ^+^-i ft 



n <i<m r (A) Ilo<i<n I 6 * + ail* -Anelfc-i,*] i<^< m l<;< m Ilo<J<n I s * + a il a ' 

A similar argument gives the type II analogue. 
Theorem 7.4. Let oq, a\, bo, ■ ■ ■ ,b m be points on the unit circle with 

1 < ao < a\ < —1, 

let lu be a point in the upper half-plane, and let , of , , . . . , (3^, t be parameters such that 

»(aj),»(a^),»(r)>0, 2(n - l)r + a + ai = ^ 

0<r<m 

where a r = a+ + /3 r = /?+ + /3~. Writing q — e(voj), we have 



lim n 



r e (g( I + 1 ) T , qir+ao i q ir+ ai ( gir+a^+af aoai; gir+aj+aj Oo / ai) giT+a +af ai /a , g 4T+Q o + a i /o Oi;p, g) 



0<i<n 

'41 (<7 a ° + a 0) /a , g^o lf g^/o 1( . . . .p 1 IKrf'W^ K,. . .rf , P ,q) 



\6(aoat 1 ]p)\ n - n ( n -V T - nao - nai ]J 



T(2(n - 1)t + a + Qi - ir)r(r) 
n . r((z + l)r)r( l r + ao)r(*r + a 1 ) 



oo.oij i<i<j<„ l<i<n 110<j"<m 
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Corollary 7.5. \19ij For any real numbers ao, <zi, 6 with —b < ao < asi, and parameters ao, cki, r w«£/i positive 
real part, 

n \ , , 11 l Xi X l\ 11 |5 + a ..|2(n-l)r+ao+ai Xl 

|ao - Oi| n( " -1)r+nao+nai-n tt r((i+ l)T)r(iT + a )r(ir + ai) 

~ |a + 6| n ( n - 1 ) T+Tlai |ai + &l n( ™~ 1)T+ ™ Qo ^ T(2(n - l)r + a + ai - ir)r(r) 

Remark. In fact (as observed in P33), the constraint that 5R(r) > is too strict, as can be seen from the fact 
that the right-hand side remains finite and positive as long as 

9ft(r) > -1/n, -5R(a )/(n - 1), -K(ai)/(n - 1). 
One can presumably weaken the conditions of Theorem 17.41 correspondingly. 

Corollary 7.6. For any real numbers ao, a\, bo, b\ with —b\ < —bo < ao < a\, and parameters ao, a\, f3o, 
(3\, t with positive real part such that ao + ct\ = /3q + /3i, we have the transformation 

where 

tt |«o + bo\ a ^\ao + frip^ 1 |oi + b \ a ^ | fll + b^~^ T(ir + ao)T{iT + aQ 
o<i * ' a ° ~ ai| 1 ~ (n ~ 1)r ~ ao ~ ai |&o - ftil-i+^-iJ-r+^o+^x T(ir + f3 )T(iT + ft) ' 

The case is similar. 

Theorem 7.7. Let ao, ■ ■ ■ , a n , bo, ... , b m , Z be points on the unit circle with 

ao < ■ ■ ■ < a n < a n +i ■= ao, 

let uj be a point in the upper half-plane, and let a^ , ■ ■ ■ , a^ , > • • • > Pm be parameters such that 

5R(a ± ),...,5R(a±) > 0, £ a r - £ ft, 

0<r<n < r < rn 

where a r — a+ + a~ , f3 r = + f3~ . Writing q — e(vuj) 7 we have 

T e {q^^;p,q) 



lim = — t 

Uo<r,s<n r e (^+ tt = a s /a r ;p, q)T e (q^ «r Z /a ■ ■ ■ a n ;p, q)T e (q^ a • • • a n /Z;p, q) 

I<Z\Z\ q"r /ar , ....... y/y/ 2 -^& r , q^a r , ....... ,^^1/2-/5+ /b r , ...;p,q) 

= JJ |g(a J /ar;p)l 1 - a '- "lg(^/ao-o rt ;p)l 1 -^°T r (Er«r) 

0<r<s<n lltKr^n 1 l"rJ 

(27r(p; P ) 2 )" /• IIo<r<nl g (^/°r;p)| a '- 1 tt gbj 

(» + ^ ^If-t:f oik " oil IW \0^ rZi ; P )\^ o 11 b 
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Remark. The fact that the n + m + 2 theta functions 

0<i<n 0<i<n 

for fixed Ilo<i<n z i s P an an n + 1-dimcnsional space is presumably relevant to finding an appropriate change 
of variables to eliminate the theta functions from the integrand. Clearly, though, the resulting computations 
would not give a trivial derivation of the above integral from a more traditional multivariate beta integral. 
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